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Abstract 

The muon is playing a unique role in sub-atomic physics. Studies of muon 
decay both determine the overall strength and establish the chiral structure 
of weak interactions, as well as setting extraordinary limits on charged-lepton- 
fiavor-violating processes. Measurements of the muon’s anomalous magnetic 
moment offer singular sensitivity to the completeness of the standard model 
and the predictions of many speculative theories. Spectroscopy of muonium 
and muonic atoms gives unmatched determinations of fundamental quantities 
including the magnetic moment ratio [igt lepton mass ratio m^lrrie^ and 
proton charge radius Vp. Also, muon capture experiments are exploring elusive 
features of weak interactions involving nucleons and nuclei. 

We will review the experimental landscape of contemporary high-precision 
and high-sensitivity experiments with muons. One focus is the novel methods 
and ingenious techniques that achieve such precision and sensitivity in recent, 
present, and planned experiments. Another focus is the uncommonly broad 
and topical range of questions in atomic, nuclear and particle physics that such 
experiments explore. 
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1. Introduction 

The muon is not a building block of ordinary matter. It’s much heavier than 
the electron but much lighter than the proton. It interacts through its electric 
charge and magnetic moment and its weak charged and neutral currents—but 
not the strong force. Positive muons can form hydrogen-like atoms with elec¬ 
trons while negative muons can form hydrogen-like atoms with nuclei. The 
muon is unstable, but sufficiently long-lived to precisely study its properties 
and sufficiently short-lived to precisely study its decays. By a quirk of nature— 
parity non-conservation—muons are produced fully polarized and when they 
decay they are self-analyzing. 

Since its discovery the muon has played a rather unique and versatile role in 
physics. In this review we discuss recent, current and near-future efforts involv¬ 
ing precision measurements of properties and decays of free muons and muonic 
atoms. The physics topics—which range from fundamental constants and basic 
symmetries, to weak nucleonic and nuclear interactions, and standard model 
tests and new physics searches—are quite diverse. These unique experiments 
are generally designed to do one thing and do it well. We aim to provide an 
experimentalist’s perspective into how these measurements are performed and 
their physics impact. 

This article is organized as follows. Section O introduces the historical 
evolution of muon experiments in sub-atomic physics. Some general comments 
on muon facilities, experimental technologies and blind-analysis procedures are 
made in Sec. ll.2[ Section [2] discusses the measurement of the muon lifetime and 
determination of the Fermi Constant. The measurements of the muon decay 
parameters and tests of the V-A theory are described in Sec. Searches for 
charged lepton flavor violating muon decays are described in Sec. Section 
covers the measurement of the muon’s magnetic and electric dipole moments, 
and their sensitivities to new particles and unknown forces. Lastly, a number 
of precision measurements involving muonium atoms and muonic {lJi~ Z) 

atoms, which determine fundamental constants and elementary interactions, are 
discussed in Secs. |§[3and|8l 

1.1. Historical overview of muon experiments in sub-atomic physics 

Table [2 summarizes the important properties and decay modes of the muon. 
The discovery of muons—new particles about 200 times more massive than 
electrons—as cosmic-ray constituents was famously made by Anderson and Ned- 
dermeyer [1] at Caltech in 1936 and Street and Stevenson [2] at Harvard in 1937. 
Their work was actually the culmination of numerous experiments conducted 
over many years, which eventually demonstrated that the highly-penetrating 
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Table 1: Summary of Measured Muon Properties and Selected Decay Rates and Limits 


Property 

Symbol 

Value 

Precision 

Ref. 

Mass 


105.658 3715(35) MeV 

34 ppb 

m 

Mean Lifetime 


2.1969811(22) x 10“® s 

1.0 ppm 

HI 

Anom. Mag. Moment 


116 592 091(63) x 10““ 

0.54 ppm 

HE] 

Elec. Dipole Moment 

dfx 

< 1.9 X 10“'^®e-cm 

95% C.L. 

la 

Branching Ratios 

PDG average 

B.R. Limits 

90% C.L. 

Ref. 


« 100% 

^ e“7 

5.7 X 10"^® 

m 


1.4(4)% 

1 

1 

+ 

1 

1.0 X 10“^^ 

HI 

^ e~h'ei'^e^e~ 

3.4(4) X lO"'^ 

e conversion 

7 X 10“^^ 

m 


component of cosmic radiation was neither electrons nor protons Q Ultimately, 
the discovery became our first evidence for generations of elementary particles 
and the hierarchical structure of the standard model. 

The first cloud chamber photograph of the decay of a muon was taken in 
1940 m- The earliest determination of the muon lifetime was made by inter¬ 
preting the “anomalous absorption” of cosmic-ray muons with decreasing alti¬ 
tude as muon disintegration [12] . The direct measurement of the muon lifetime 
was made shortly afterwards by recording the time intervals between stopping 
muons and decay electrons using cosmic rays [13] . These early measurements— 
involving in-flight and stopped muons—afforded a decisive test of time dilation 
for moving particles. Nowadays the precision measurement of the muon life¬ 
time provides the best determination of the Fermi constant Gi?, the quantity 
governing the universal strength of weak interactions. 

The /i ^ euu decay scheme was established through cosmic-ray measure¬ 
ments of the decay electron energy spectrum. It was demonstrated—via the 
continuum distribution and the energy endpoint—that muons decay into three 
or more particles with small or zero masses m- Fermi’s theory of nuclear 
/3-decay was consequently expanded to accommodate fi evv decay thus af¬ 
fording the first glimpse of weak universality. These early experiments—along 
with theoretical work on the tensor structure of the current-current interaction 
[15] [16] —were the beginnings of our modern precision studies of the muon decay 
parameters as a valued probe of the weak force. 

In a profound paper in 1956, Lee and Yang suggested that the discrete 
symmetry of parity might be violated in the weak interaction HZ]. One impor¬ 
tant prediction of parity non-conservation concerned the by-then well-known 
TT ^ /i ^ e weak decay chain. First, the non-conservation of parity in ir ^ jau 
decay would cause muons to be polarized along the muon momentum axis. Sec¬ 
ond, the non-conservation of parity in /i ^ decay would cause electrons to 
be emitted anisotropically about the muon polarization axis. Soon afterwards 


^The existence of highly-penetrating cosmic rays was known since the work of Bothe and 
Kolhorster in 1929. In 1933, Kunze [3] noted a “particle of uncertain nature” in his investi¬ 
gations of cosmic rays using ionization chambers. 
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the non-conservation of parity in tt ^ /i ^ e decay was reported by Garwin 
et al. ng, and Friedman and Telegdi m, following their observation of large 
decay-electron anisotropies. Parity non-conservation is now codified in the V-A 
structure of the weak currents [20l |2T] which imparts a distinctive angle-energy 
correlation on the decay electrons with the muon polarization. Since this early 
work, the precision measurement of angle-energy correlations has enabled in¬ 
creasingly precise tests of V-A theory. 

The hypothesis of a weak interaction mediated by a force-carrying boson 
was originally introduced in 1940 [22]. The early theories involving weak bosons 
predicted the existence of /i ^ ey decay at the level of about 10“^ (the mecha¬ 
nism for /i ^ ey involved the neutrino emitted by the muon being absorbed by 
the electron). By the late 1950s the /i ^ ey experimental limit was orders-of- 
magnitude below this theoretical prediction. The crisis concerning /i ^ ey decay 
led Pontecorvo [23] and others to postulate the existence of two distinct neu¬ 
trino types; an electron-flavored neutrino and a muon-flavored neutrino. Since 
the genesis of lepton flavor in /i ^ ey, increasingly delicate searches for rare 
processes including /i ^ ey and ji ^ e conversion have continued to shape our 
understanding of flavor. 

The aforementioned Garwin et al. experiment on parity non-conservation 
in muon decay also yielded the first measurement of the magnetic moment 
of the positive muon. Their result for the gyromagnetic ratio g = +2.00 ± 
0.10 was a demonstration that the muon was a structureless, spin-1/2 Dirac 
particle. Subsequently—through increasingly sophisticated measurements that 
utilize the possibilities of polarization and polarimetry of muons in tt ^ /i ^ e 
decay—the determination of the anomalous part {g — 2) of the muon magnetic 
moment has been measured to an astonishing sub-part-per-million level. The 
anomalous moment arises through quantum vacuum fluctuations that accrue 
from all particles of nature—both known and unknown—and thereby affords a 
unique window on new physics at high-energy scales 

The simplest atom involving muons is muonium, a pure QED bound state 
of sometimes dubbed the “perfect atom.” Because both constituents are 

point-like leptons the muonium energy levels are completely free of perturbations 
arising from nuclear size effects. Its formation was first identified by detecting 
the characteristic Larmor precession frequency of polarized muonium formed 
when muons stop in certain gases. Following its discovery the measurement of 
the hyperfine splitting of the muonium ground state—by inducing microwave 
transitions between hyperfine states in external magnetic fields—was developed 
by Hughes and co-workers (for details see [24]). This work has provided the most 
precise determinations of the muon-to-proton magnetic moment ratio and 

the muon-to-electron mass ratio m^lrrie as well as important tests of quantum 
electrodynamics. 

When negative muons are brought to rest in matter they undergo atomic cap¬ 
ture and form muonic atoms. These atoms are hydrogen-like systems where— 
unimpeded by the exclusion principle—the muon cascades from an initially high 
principal quantum number state to the 15' atomic ground state. The existence of 
such atoms was first discussed in Refs. [25] [26] that demonstrated the timescale 
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for formation of muonic atoms was much shorter than the muon lifetime. Com¬ 
pared to ordinary atoms, the muonic atom radii are (m^/mg) times smaller and 
the energy levels are (m^/mg) times greater. Consequently, the overlap between 
the muon orbits and the nucleus is much larger than in ordinary atoms and the 
energy levels can be significantly perturbed by the nuclear charge distribution. 
Precision spectroscopy of muonic atoms thus became a workhorse for studies of 
nuclear charge radii and electromagnetic moments. In the simplest case of a /ip 
atom, the recent measurements of the tiny energy splitting (Lamb shift) between 
the muonic 2 Si /2 — 2 Pij 2 orbitals have provided the most precise determination 
of the proton charge radius. 

After muonic atoms are formed they disintegrate by either muon decay or 
nuclear capture /i“[A, Z] ^ [A, Z — l]u. The first experimental evidence for 
muon capture was observations of differing electron yields when stopping neg¬ 
ative and positive muons in carbon and iron EH- This work demonstrated 
the weak nature of muon capture on atomic nuclei thus dispelling the notion 
that cosmic-ray muons were the force-carriers of the strong interaction, and 
ultimately leading to early ideas of weak universality in muon capture, muon 
decay and beta decay. Modern measurements of muon capture are investigating 
weak nucleonic and nuclear interactions that address subjects which range from 
standard model symmetries to fundamental astrophysical processes. 

1.2. Common features of preeision muon experiments 

Muon beams are derived from pion decays, the pions being produced in the 
nuclear collisions between an accelerated beam and a fixed target. There are 
three basic types of muon beam lines involving so-called surface, cloud and de¬ 
cay muons. In surface pA beams the muons originate from at-rest decay of 
7r+ stops in the surface layer of the production target. The resulting muons 
are mono-energetic 29.8 MeV/c, 100% longitudinally polarized, and because of 
the localized source have a sharp focusIn cloud beams the muons originate 
from in-fiight decays of parent pions in the region between production target 
and the first bending magnet of the secondary beam line. The resulting muons 
arise from “forward-decays” and “backward-decays” of pions and therefore the 
resulting polarization is considerably lower than surface beams. In decay beams 
an upstream section of the beam line selects the parent pion momentum and 
a downstream section of the beam line selects the daughter muon momentum. 
The resulting muons are typically highly polarized and free of electron contam¬ 
ination. 

Intense muon sources at the Paul Scherrer Institute (PSI) in Switzerland 
and TRIUMF in Canada are examples of muon beams based on high current, 
medium energy, proton cyclotrons that provide an essentially continuous beam 
(with the micro-time structure of the cyclotron radio-frequency). Facilities at 
Fermilab in the U.S. and J-PARC in Japan are examples of muon beams based 


^Only surface beams are available as tt stops form pionic atoms and rapidly undergo 
nuclear capture. 
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on lower current, higher energy, proton synchrotrons and provide a pulsed beam 
with typical repetition rates of tens of Hertz. Additionally, novel muon sources 
for ultra-cold muons, based on formation and ionization of muonium atoms, 
and ultra-intense muons, based on capture and transport by superconducting 


solenoids, are under development for future muon experiments (see Secs. 5.3.2 
and 4.2.3| respectively). 

Today’s precision muon experiments—which require both enormous statis¬ 
tics and extraordinary limits on possible biases from systematic effects—are 
benefiting from advances in areas including radiation detectors, readout elec¬ 
tronics, computer hardware and software infrastructure. 

Challenges for detector design include stringent requirements on timing, 
tracking and calorimetry in high-rate environments. Ultra-low mass tracking 
chambers and modern silicon detector technology have been used or are being 
developed for high-rate, high-precision tracking applications. High-density elec¬ 
tromagnetic calorimeters (liquid xenon, lead fluoride, lead tungstate) have been 
built or are being constructed for both good energy resolution and fast-timing 
applications. Modern silicon photomultipliers are enabling readout of scintil¬ 
lation and Cherenkov light in magnetic fields and restricted geometries. And 
sophisticated NMR-based field measurements and finite-element field modeling 
are being utilized for precision magnetometry. 

Commercial applications for real-time measurement, processing and distri¬ 
bution of massive data-streams have driven hardware development including 
high sampling-rate waveform digitizers, field programmable gate arrays and 
graphical processing units. The digital capture of detector signals offers in¬ 
valuable opportunities for systematics investigations in precision measurements 
and waveform digitizers have been used or are being deployed in many recent, 
present and future experiments. Field programmable gate arrays— i.e. high- 
performance user-customizable integrated circuits—are now common for special 
purpose tasks in timing, trigger and control logic. A hybrid architecture of mul¬ 
ticore CPUs and teraflop-performance GPUs is being developed for real-time 
readout and processing in the Fermilab muon {g — 2) experiment. 

Recent precision muon experiments have stored datasets of hundreds of ter¬ 
abytes of raw data and future experiments will store datasets of many petabytes 
of raw data. This scale of data analysis and data simulations has required grid 
computing facilities such as WesGrid (Canada) and the National Center for 
Supercomputing Applications (U.S.). The extreme demands on high statistics 
and understanding systematics are pushing applications of GEANT (GEometry 
ANd Tracking software toolkit) into new territories. 

Increasingly, modern measurements are carried out using so-called “blind- 
analysis” techniques. A simple example for measurements that use precision 
oscillators for timing is to prescriptively de-tune the oscillator by a small off¬ 
set from its nominal setting during the data taking; the offset is unknown to 
anyone analyzing the data. The data analysis time unit is then a somewhat ar¬ 
bitrarily defined “clock tick.” When the analyses are complete, the clock ticks 
are converted into physical time units, and the unblinded result is revealed. In 
efforts where rare events are investigated, blinded regions must be established 
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Figure 1: a). Tree-level Feynman diagram for ordinary muon decay in Fermi’s current-current 
interaction, b). Tree-level standard model diagram for ordinary muon decay indicating the 
W-boson mediated weak interaction between the leptonic currents. 


in which the rare events would be found. Backgrounds and cuts are determined 
from allowed events outside and near the forbidden box. Once the analysis is 
complete, the blinded box is opened to see if any events have survived. These 
disciplined procedures provide a needed level of integrity to the experiments. 


2. Muon lifetime 


2.1. Fermi constant Gp 

The strength of the weak interaction is governed by the Fermi constant 
Gp- The roots of are Fermi’s theory—based on an analogy between the 
emission of an electron-neutrino pair by a radioactive nucleus and a photon 
by a charged particle—of a current-current weak interaction. Of course, since 
1934, our modern understanding of weak interactions has evolved to incorpo¬ 
rate parity-violating V — A currents and the massive W and Z gauge bosons. 
However, the constant Gp and Fermi interaction have survived as a convenient, 
low energy, effective theory of the weak sector in the standard model (and pre¬ 
sumably any successor). 

Within the standard model the Fermi constant (see Fig. is given by 


Gp 

7f 


8M^ 


1 + ^5 


( 1 ) 


where l/M^ represents the tree-level propagator corresponding to VF-boson 
exchange and g the weak coupling. The term Vi incorporates the higher- 
order electroweak interaction corrections [28]. The factors of \/2 and 8 in Eqn. 

are reminders of the origins of the Fermi constant in a vector current - vector 
current weak interaction. 

By far the best determination of the Fermi constant is obtained by the 
measurement of the positive muon lifetime, r^. Experimentally, intense beams 
of low-energy muons are nowadays available and the 2.2 /rs muon lifetime with 
its associated decay electrons are nicely suited to precision measurements of 
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time distributions. Theoretically, because muon decays—/i ^ ez/P, /i ^ ez/z /7 
and li eeevv —are pure leptonic weak interactions, their interpretation is 

unambiguous. 

The determination of the Fermi constant Gf from the muon lifetime 
represents a reference point for subatomic physics. It permits the testing of 
weak universality, for example, through precision measurements of leptonic tan- 
decays. It enables the determination of weak mass-mixing angles, for example, 
through precision measurements of neutron decay. Moreover—together with 
the fine structure constant a and the Z gauge boson mass Mz —it completely 
determines the electroweak sector of the standard model and enables searches 
for new forces and particles. 

2.2. Experimental approaehes to measuring 

As already discussed, after the discovery of the muon by Anderson and 
Neddermeyer [T] and Street and Stevenson [ 2 ] the earliest measurements of the 
lifetime—for stopped muons and in-flight muons—were important in verifying 
the time dilation of moving particles. 

By the beginning of this century the Particle Data Group world average 
of the muon lifetime was = 2.19703 ± 0.00004 /is or 18 parts-per-million 
(ppm) [29]. The world average was largely determined by three measurements: 
Giovanetti et al. EO], Bardin et al. EH and Balandin et al. [32], that were con¬ 
ducted in the seventies and the eighties. The experiments of Giovanetti et al. 
and Balandin et al. used low-rate continuous beams in order to insure the arrival 
and decay of muons occur one-by-one, i.e. avoiding any incorrect assignment 
of daughter electrons with parent muons. At higher rates, if only the previous 
stop was associated with a particular electron the measured lifetime would be 
distorted, or, if all the neighboring stops were associated with a particular elec¬ 
tron a random background would be incurred. Such one-by-one measurements 
therefore limit the collection of decays to roughly 10 ^^ and the uncertainty on 

to roughly 10 ppm. 

The two most recent measurements of the positive muon lifetime—the FAST 
experiment [33] and MuLan experiment [5]— were specifically designed to cir¬ 
cumvent the statistical limitations of one-by-one measurements. The FAST 
approach involved an active pixelated target in order to reconstruct muon- 
electron vertices and thereby correctly associate each decay electron with its 
parent muon. In principle therefore, multiple muons could be simultaneously 
stopped and multiple decays could be simultaneously recorded, without los¬ 
ing the parent-daughter association. The MuLan approach involved a time- 
structured muon beam in order to first prepare a “radioactive source” of muons 
and afterwards measure the “emanating radiation” of electrons. In this scheme 
no association of a particular daughter electron with a particular parent muon 
is necessary; the observed lifetime of a radioactive source is independent of the 
source preparation. 
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Figure 2: Cutaway diagram of the MuLan experiment indicating an incoming muon and 
outgoing positron and showing the beam pipe (1), stopping target (2), optional magnet array 
(3), scintillator detector (4) and beam monitor (5). Figure courtesy MuLan collaboration. 


2.2.1. MuLan experiment 

The MuLan experiment [5] was conducted at PSI. The setup—which com¬ 
prised an in-vacuum stopping target and fast-timing finely-segmented positron 
detector—is shown in Fig. The setup was designed to both accumulate the 
necessary quantity of muon decays and minimize distortions arising from muon 
spin rotation in the stopping target and electron pile-up in the detector array. 

To reach the necessary statistics of 10^^ decays the experiment relied on 
a time-structured surface-muon beam (see Fig. |^. The experiment involved 
cycles of 5 /rs duration beam-on periods to accumulate stopped muons and 
22 /is duration beam-off periods to measure decay positrons. The arrangement 
permitted an average rate of stopped muons of 1-2 x 10^ s“^; much more than 
permissible in a one-by-one lifetime measurement. 

Surface muon beams are nearly 100% longitudinally polarized. The positrons 
emitted in muon decay are distributed asymmetrically about the /i-spin axis with 
high-energy positrons preferentially emitted in the spin direction and low energy 
positrons preferentially emitted opposite the spin direction. The spin vectors of 
stopped muons both process and relax in the local magnetic field of the target 
material; a phenomenon known as muon spin rotation. /iSR yields a time- 
dependent muon-ensemble polarization and thereby a time-dependent decay- 
positron angular distribution. It results—when detecting positrons in specific 
directions—in a geometry-dependent modulation of the exponential decay curve 
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Figure 3: The 5-yLiS beam-on, 22-/i,s beam-off, time structure of the MuLan experiment. The 
upper panel shows the muon arrival time distribution and the lower panel shows the decay 
electron time distribution. 


by the /iSR signal. 

Two combinations of stopping target material and environmental magnetic 
fields were used. One configuration involved a magnetized ferromagnetic foil 
with a roughly 0.5 T internal field orientated perpendicular to the beam axis. 
Another configuration involved a nonmagnetic quartz crystal with a roughly 
80 G external field orientated perpendicular to the beam axis. Muon stops in the 
ferromagnetic target were mostly diamagnetic /i+ ions exhibiting a precession 
frequency of 13.6 kHz per Gauss while muon stops in the quartz target were 
mostly paramagnetic /i+e“ atoms exhibiting a precession frequency of 1.39 MHz 
per Gauss. The targets were mounted in the beam vacuum to avoid either 
multiple scattering or stopping in upstream detectors or vacuum windows. 

The 5 /is duration of muon accumulation was important in dephasing the 
spins of the polarized stopping muons in the transverse magnetic fields of the 
two target configurations. In both cases the precession frequencies of either 
/i+ ions or atoms was sufficient to reduce the transverse polarization by 

roughly a thousand-fold. The residual polarization was ultimately limited by 
the alignment accuracy between the transverse field and the beam polarization, 
i.e. the presence of a small longitudinal component of the muon polarization. 

The positron detector comprised a nearly 47r array of 170 fast-timing, double¬ 
layered, scintillator detectors arranged in a truncated icosahedron (soccer ball) 
geometry. The high granularity and fast timing characteristics were important in 
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minimizing the incidence of positron pile-up. The forward-backward symmetry 
of the Stt solid angle detector also suppressed the imprint of /iSR on the time 
distribution of the decay positrons]^ Therefore the MuLan /iSR signal is strongly 
suppressed and largely determined by the the non-uniformity of the detection 
efficiency. 

Analog signals from detector elements were digitized using fast-sampling 
ADCs and digitized “islands” of contiguous samples of above-threshold signals 
were identified by FPGAs. A distributed data acquisition enabled the storage 
of all above-threshold signals from the positron detector. 

A temperature-stabilized crystal oscillator was used as the timebase for the 
fast-sampling ADCs. The collaboration was blinded to the exact frequency of 
the timebase during the data taking and the subsequent analysis. Only after 
completing the entire analysis was the frequency unblinded and the lifetime 
revealed. 

In analyzing the data the digitized islands were first fit to pulse templates to 
determine the times and energies of individual pulses. A software-defined min¬ 
imum amplitude was applied to distinguish the minimally ionizing positrons 
from low-energy backgrounds and a software-defined minimum deadtime was 
applied to establish an explicit resolving time between neighboring pulses. The 
time distribution of coincident hits between inner-outer tile pairs was then con¬ 
structed. 

After the application of small data-driven corrections for positron pileup and 
gain changes, the time distributions of coincident hits were fit to extract r^. The 
ferromagnetic target data showed no evidence of /rSR effects and was fit to the 
function 

N{t) = + B 

where the time-independent background B arose from both cosmic rays and the 
imperfect beam extinction during the measurement period. The quartz target 
data showed a clear /rSR signal and was therefore fit to a modified function that 
incorporated both longitudinal- and transverse-field /iSR effects. The fits were 
performed for different software deadtime and extrapolated to zero software 
deadtime to obtain r^. 

The values of obtained from the two targets were in good agreement and 
yielded a combined result of 

= 2 196 980.3 ± 2.1(stat) ± 0.7(syst) ps; 

an overall uncertainty of 2.2 ps (1.0 ppm) and thirty-fold improvement over 
earlier generations of experiments. 

2.2.2. FAST experiment 

The FAST experiment [33] was also conducted at PSI. The setup—which 
consisted of beam defining counters and a pixelated stopping target—is shown 


detector with perfectly forward-backward symmetry about the muon spin direction 
would display no /rSR signal in the positron time distribution. 
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Figure 4: Horizontal cross section of the FAST experiment showing (from left to right) 
the wedge-shaped degrader, beam defining scintillators and pixelated stopping target. The 
pattern of hits from a tt —)■ /r —)■ e decay chain are also indicated. Figure courtesy FAST 
collaboration. 


in Fig[^ 

The approaches of FAST and MuLan to measuring were very different. In 
FAST the segmented target was used to associate the daughter positrons with 
parent muons by identifying the coordinates of the muon stop and the trajectory 
of the decay electron. 

Additionally, the FAST experiment employed a stopping pion beam. When 
pions were stopped the subsequent muons from tt ^ jiu decay and subsequent 
positrons from /i ^ euu decay were detected, thereby identifying the entire 
TT ^ ^ e decay chain. Since the decay muons from stopped pions are emit¬ 

ted isotropically, the overall ensemble of stopped muons had essentially zero 
polarization. This minimizes the /iSR distortions of the decay-positron time 
distribution. 

To widely distribute the stopping pions over the segmented target the beam 
was defocused at the target location and directed through a wedge-shaped de¬ 
grader. The degrader varied the penetration depth of pions according to their 
vertical coordinate. Beam counters were situated both upstream and down¬ 
stream and provided the identification of stopping pions and through going 
electrons. A permanent magnet comprising two planes of ferrite blocks pro¬ 
vided a 80 G transverse magnetic field. 

The segmented stopping target comprised a 32 x 48 pixelated array of 4 x 
4 X 200 mm^, vertically orientated, plastic scintillator bars. The light from in¬ 
dividual bars was readout by multi-anode photomultipliers. After amplification 
and discrimination the hit times of individual pixels were recorded by multihit 
TDCs. The setup employed two discriminator thresholds, a higher-level thresh- 
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old capable of identifying the heavily ionizing pions / muons that was used in 
the trigger logic and a lower-level threshold capable of identifying the minimally 
ionizing positrons that was used in the data analysis. 

In normal data-taking the TDC information was readout after the identifi¬ 
cation of a pion stop and a corresponding decay muon in the pixelated target. 
A level-one trigger distinguished the incident pions from other beam particles. 
A level-two trigger identified the sequence involving the prompt pion signal and 
the delayed muon signal (it used an FPGA array to perform the trigger decision 
within several microseconds). The range of the 4.1 MeV muons from the i: ^ fiv 
decay was approximately 1.5 mm and therefore the muon stop was located in 
either the same pixel or an adjacent pixel to the pion stop. On fulfilling the 
above trigger the pion stop pixel was used to define a 7 x 7 pixel region and a 
—8 to +22 /is time window for the selective readout of the TDC modules. 

Data were collected at incident beam rates of roughly 160 kHz and stopping 
pion rates of roughly 80 kHz, and yielded a sample of approximately 1.0 x 10^^ 
fi evv decay events. To handle the large data volume the TDC information 
was both readout and also processed in real-time. The real-time processing used 
lookup tables to identify the characteristic topologies of pixel hits corresponding 
to a positron emanating from the muon stop (the topologies allowed for pixel 
inefficiencies). 

The identification of tt ^ /i ^ e events yielded the corresponding times C, 
and te of the pion, muon and positron and thereby the time differences between 
the stopped pion and the decay electron, — C, and between the stopped 
muon and the decay electron, This real-time processing accumulated 

both global te — and te — histograms for all tt ^ /i ^ e events as well 
as sub-sets for different pion coordinates, positron topologies, etc. Both the 
measured tg — C and tg — distributions exhibit the muon lifetimej^ 

The time distributions showed the exponential decay curve as well as time- 
independent and cyclotron-RF correlated backgrounds. The two backgrounds 
arose from random backgrounds and beam particles that generated fake positron 
topologies and thereby tt ^ /i ^ e events. The time distributions showed no 
evidence of /iSR effects. 

The RF-correlated background was handled by rebinning data with the cy¬ 
clotron RF period. After the rebinning procedure a maximum likelihood fit to 
the decay curve gave the final result 

= 2 197 083 ± 32(stat) ± 15(syst) ps 
with an overall uncertainty of 35 ps (16 ppm). 

2.3. Negative muon lifetime 

The above experiments involved measurements of the lifetime of the pos¬ 
itively charged muon. As discussed elsewhere, for stopped negative muons, in 


"^The te — + distribution is a single decay curve with lifetime T/j, for time intervals much 
greater than Ttt- 
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addition to muon decay the process of muon capture is possible. Measurement 
of the negative muon lifetime in a particular stopping material therefore deter¬ 
mines a total disappearance rate, i.e. the sum of the decay rate and the capture 
rate. The contribution of the capture rate to the disappearance rate varies from 
roughly 0.1% in hydrogen isotopes to greater than 90% in heavy elements. 

In principle the lifetime of free negative muons can be determined by the 
combined measurements of the disappearance rate and the capture rate for a 
specific stopping target. For muon stops in hydrogen and deuterium there exist 
both measurements of the muon disappearance rates {^^D) via the decay electron 
time spectra and the total muon capture rates (A^) via the capture neutron 
yields. Thus the combination of these measurements allows the determination 
of the free negative muon lifetime using 

1/^/x- = - Ac- (2) 

Unfortunately, the procedure is complicated by muon atomic and molecular 
processes that occur in isotopes of hydrogen. Consequently, to directly extract 
the lifetime r^- the values of Ad and Ac must be determined for the same 
muonic atomic and molecular state populations. We note that together—the 
recent measurement of the muon disappearance rate [34] and the earlier mea¬ 
surements of the muon capture rate [35l [36] from the singlet fi~p atom—yield 
using Eqn. a negative lifetime that is entirely consistent with the positive 
muon lifetimel£] 

2.4 ‘ Results for the muon lifetime 

The uncertainty of the MuLan measurement of the muon lifetime is about 
30-40 times smaller than the previous generation of lifetime experiments and 
about ten times smaller than the published result from the FAST experiment. 
While the MuLan result is in reasonable agreement with the earlier experiments 
there is some tension between the MuLan experiment and the FAST experiment. 
The weighted average of all results gives a lifetime = 2 196 981.1 ±2.2 ps with 
a chi-squared value that is dominated by the 2 a difference between MuLan and 
FAST. 

As discussed in Sec. the accurate knowledge of is important to preci¬ 
sion measurements in muon capture. The MuCap experiment |37| has recently 
measured the /i“p singlet capture rate A^ and the MuSun experiment [38] is 
currently measuring the jii~d doublet capture rate A^^. Both experiments derive 
these capture rates from the tiny difference (A^ — A) between the positive muon 
decay rate (Aq = 1/t^+) and the muonic atom disappearance rates (A = 1/r^z)- 
Because the capture rates are very small, very precise determinations of both the 
muonic atom disappearance rates and the free muon decay rate are necessary. 


^Small corrections for ii~p bound state effects and ppfi molecular formation are required 
to extract of the negative muon decay rate from the measured values of the disappearance 
rate and the capture rate. See Ref. m for details. 
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Prior to this work the muon decay rate was known to 9 ppm or 5 s“^, 
an uncertainty that limited the capture rate determinations from disappearance 
rate experiments. The recent work on the muon lifetime—yielding a precision 
of 1 ppm or 0.5 s~^ in —has eliminated this source of uncertainty. 

2.5. Results for the Fermi eonstant 

As mentioned, the value of the Fermi constant Gf can be extracted from the 
muon lifetime r^. Until recently the uncertainty in the theoretical relationship 
between the Fermi constant and the muon lifetime was about 15 ppm. 

Recent work by van Ritbergen and Stuart [39l |40l [41] and Pak and Czar- 
necki [42] have reduced this theoretical uncertainty to about 0.14 ppm. Van 
Ritbergen and Stuart were first to calculate the muon lifetime with 2-loop order 
QED corrections. Pak and Czarnecki were subsequently responsible for treating 
corrections from the mass of the electron. 

The van Ritbergen and Stuart relation between the Fermi constant Gf and 
the muon lifetime —derived using the V-A current-current Fermi interaction 
with QED corrections evaluated to 2-loop order—yields 


Gf 


I192f^ 


1 -h -1- + Ag(^) 


(3) 


where is the measured muon lifetime, is the measured muon mass, and 
Ag^^^ and Ag^^^ are theoretical corrections with Ag^^^ accounting for 
the effects of the non-zero electron mass on the muon-decay phase space and 
Ag^^/^^ accounting for the contributions of the l-/2-loop radiative corrections 
to the decay amplitude. 

Using Eqn. with the MuLan value for the muon lifetime [5], the 2010 
CODATA recommended value of the muon mass [4], and the theoretical 
corrections Ag^^\ Ag^^^ and Ag^^^ gives 


Gf = 1.166 378 7(6) x 10“^ GeV“^ (0.5 ppm). 


The result is a thirty-fold improvement over the Particle Data Group value prior 
to the recent theoretical work and the MuLan / EAST experiments. 

The electroweak sector of the standard model involves three parameters, the 
two gauge coupling constants g, gf and the Higgs vacuum expectation value v. 
Their values are fixed by measurements of the fine structure constant, a, Eermi 
coupling constant, Gf, and Z boson mass, M^. Gonsequently, the thirty-fold 
improvement in the determination of Gf, together with other improvements in 
the determinations of a and M^, now permit improved tests of the electroweak 
sector of the standard model. 

Within the standard model the quantities a, Mz and Gf are related to 
other fundamental quantities that include the charged weak boson mass Mw 


®This value for Gp is also listed in the 2014 compilation of the physical constants by the 
Particle Data Group m- 
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and weak mixing angle Ow- Such SM relationships have radiative corrections 
that impart sensitivities to the top quark mass rrit and the Higgs boson mass 
rrih- Historically, precision electroweak data was important in constraining the 
then-unknown masses of the top quark and the Higgs boson. 

The top quark was discovered in 1995 at Fermilab by the CDF and DO 
experiments mma and the Higgs boson was discovered in 2012 at CERN by 
the ATLAS and CMS experiments [46l [47] . The spectacular agreement between 
the forecast and the measurement of top mass rrit = 173 GeV/c^ was a huge 
triumph for the predictive power of the standard model. Through the SM 
relations for the Higgs sector [48] 

Gf = 
rrih = 2z/^A 

our knowledge of Gf and rrih are sufficient to determine the two parameters— 
the vacuum expectation value v and self-interaction parameter A—of the SM 
Higgs potential. 

The Fermi constant obtained from the muon lifetime is the anchor of the 
universality tests of the weak force. Comparison between the purely leptonic 
decays of the tan and the muon—z.e., r ^ r ^ euu and fi evv —are 
natural opportunities for testing the universality of leptonic weak interactions 
across the three generations. Using the available data on the r lifetime and its 
purely leptonic [ivv and evv branching ratios the leptonic universality of weak 
interactions has been demonstrated to the levels of several parts-per-thousand 
[43]. 

The Fermi constant is also the metric for extracting the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements which specify the flavor mixing in weak in¬ 
teractions of quarks. The concept of flavor mixing of quark states was originally 
introduced by Cabibbo to explain the different strengths of strangeness-changing 
and strangeness-conserving weak interactions and thereby save weak universal¬ 
ity. Modern tests of weak universality for three quark generations are therefore 
based on testing whether the sum of all couplings of an up-type quark (up, 
charm or top) to all down-type quarks (down, strange and bottom) is equal to 
unity. For the case of the up-quark, one obtains [49] 

\Vud\^ + iKaP + iKbP = 0.99978(55) 

which demonstrates weak universality between the quark-lepton sectors at the 
level of about 0.6 parts-per-thousand. 

3. Muon decay 

3.1. V-A structure of the weak interaction 

Herein we discuss the decay /i+ ^ and the investigation of the 

Lorentz structure of the leptonic weak charged current. Through combina¬ 
tion of its purely-leptonic nature—permitting rigorous theoretical analysis— 
and its not-too-short, not-too-long lifetime—permitting rigorous experimental 
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Figure 5: Left: Relative probability in arbitrary units of the emitted positron in ordinary 
muon decay as a function of the muon spin direction and the emitted positron reduced 
energy y for the V-A weak interaction. Right: Projections of the relative number, N(y), and 
asymmetry, A(y), vs. positron reduced energy. 


analysis—the decay ^ is the quintessential test of the V-A structure 

of the weak interaction. 

The measurable observables of emitted positrons in /i+ ^ e+z^P decay are the 
positron energy and directional distributions, shown in Fig. and the positron 
polarization. The V-A interaction dictates the precise form of the energy- 
angle distribution and the polarization observables which include the parity 
non-conserving angular asymmetry and longitudinal polarization. In general 
the presence of non-V-A currents would change these energy, angle and polar¬ 
ization distributions. 

Based upon a general current-current interaction with possible scalar, pseu¬ 
doscalar, tensor, vector and axial terms the positron energy, angle and po¬ 
larization distributions may be parameterized by nine parameters (see Refs. 
[ISlIIiEo]). The differential probability for emitting a positron with a reduced 
energy y at an angle 0 from the muon spin axis is given by 

1 /- 

dydcosd ^ V2/' - 2/0 (Fisiy) + cos ePAsiv)) (4) 


where = \P^\ is the muon polarization, VFe/i = = ^2.8 MeV 

is the maximum positron energy, y = is the reduced energy, and yo = 

= 9.7 X 10“^. The y-dependent functions 

Fis{y) = 2/(1 -y) + ^p( 42/^ - 3y - yl) + 7?yo(l - y) (5) 


FAs{y) 




9 ? 1 


2, 

9+3^ 


4y - 3 + (yi - _ i) 


(6) 
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describe the isotropic and anisotropic parts of the energy spectrum. They are 
governed by the decay parameters p, p, (^, and S—conventionally termed the 
Michel parameters. For unpolarized muons, the two parameters p and p com¬ 
pletely determine the energy distribution with p governing the high-energy part 
and 77 governing the low-energy part of the spectrum^ For polarized muons, the 
two parameters ^ and 5 additionally determine the angular distribution with ^ 
governing the energy-integrated positron asymmetry and S governing the en¬ 
ergy dependence of this angular asymmetry. The standard model values of the 
Michel parameters are p = 0.75, p = 0.0, ^ = 1.0 and 6 = 0.75. 

The energy-angle dependent positron polarization vector Pe{y,0) may be 
decomposed into a longitudinal component Pl and two transverse components 
Ptii Pt 2 according to 

Pe{y,0) = Pl ■ Z + Pti ■ {z P^) z + Pt 2 ■ {z x P^) (7) 


where i and are unit vectors along the positron momentum and the muon 
polarization, respectively. The transverse polarization Pti lies in the decay 
plane of the positron momentum and the muon spin while the transverse po¬ 
larization Pt 2 lies along the perpendicular to this decay plane. Note that the 
transverse polarization Pti is a (time-reversal) T-conserving observable while 
the transverse polarization Pt 2 is a T-violating observable. 

The energy-angle dependence of the longitudinal polarization Pl is given by 


p / n\ ^ Pip{y) + Pti. COs9FAp{y) 
^ ’ Fis{y) + Pn cos OF AS {y)' 


The two additional y-dependent functions are 


(8) 


Fip{y) = ^Vy^ - yl (-2y + 2 + yi -yt) + 4^ (py-A + 

(9) 

FAp{y) = ^ |c" (2^ - y -yl) + (4^ - Sy - yl) + 2-q"{l - y)yo^ 

( 10 ) 

and involve two further decay parameters—and Neglecting the positron 
mass and radiative corrections, the standard model values of the decay param¬ 
eters = 1.0 yield an energy-angle independent longitudinal polarization 

Pl = 1. 

The T-conserving polarization Pti and its energy-angle dependence is given 


by 


PTi{y,0) 


P^ sin OFti jy) 
Fis{y) + Pn cos OPAsiy) 


( 11 ) 


yl 


^The contribution of 77 is order yo and consequently very small. 

®An additional parameter y" also enters FAp{y) but its role is highly suppressed by the 
factor yo. 
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and the T-violating polarization Pt 2 and its energy-angle dependence is given 
by 


PT2{y,0) 


P^SUlBFT2{y) 
Fisiy) + Pn cos OFAsiv )' 


( 12 ) 


The two additional y-dependent functions are 


Friiy) = Y2 




^" + I2[p- 


FT2{y) = \^y'^ - Vo 


(1 - y)yo - ‘iy{y‘^ - yl) + r]"{-‘iy‘^ + % - yl) 

(13) 

(14) 


3^(i-y)+2§yr^ 


The T-conserving polarization Pti and its energy-angle dependence are de¬ 
termined by the decay parameters 77 and rj" and the T-violating polarization 
Pt2 and its energy-angle dependence are determined by the decay parame¬ 
ters a'/A and /A. The standard model values of these parameters are all 
zero; r] = r]" = a' jA = jA = 0. They yield—when neglecting the positron 
mass and radiative corrections—an energy-angle independent transverse polar¬ 
ization Pti = Pt 2 = 0. On accounting for positron mass effects the average 
T-conserving transverse polarization component is {Pti) = —3 x 10 “^ [51]. 

The standard model and measured values of the decay parameters are sum¬ 
marized in Table Until recently our knowledge of the muon decay parameters 
was based on experiments from twenty, thirty and forty years ago. For example, 
the best determination of the energy spectrum parameter p was due to Derenzo 
et al [52] and the best determinations of the angular distribution parameters ^ 
and 6 were due to Balke et a/. [53], Beltrami et a/. [54] and Jodidio et al ESI- 

The past decade has seen a substantial improvement in the experimental 
knowledge of the decay parameters of both the /i+ ^ e^uu energy-angle distri¬ 
bution and the various polarization observables. A precision measurement of the 
positron energy-angle distribution in polarized muon decay was conducted by 
the TWIST collaboration [56] |57] at the TRIUMF cyclotron. It led to roughly a 
factor ten improvement in the determination of the decay parameters p, ^ and 6 . 
Other precision measurements of both the transverse polarization [51] and the 
longitudinal polarization [58] of the decay positrons were recently conducted at 
PSI. They have improved our knowledge of the longitudinal polarization param¬ 
eter by roughly a factor of ten and the transverse polarization parameters 
by roughly a factor of three. 


3.2. TWIST experiment 

TWIST—the TRIUMF Weak Interaction Symmetry Test experiment [56] 
l57] —was designed to measure with unprecedented accuracy the energy-angle 
distribution of decay positrons from polarized muons. The setup—including 
the muon beam, stopping target and positron spectrometer—is shown in Fig. 
[^ The experiment involved stopping a highly polarized muon beam in a very 
thin target foil and detecting the resulting decay positrons in a low mass, large 
acceptance, high resolution magnetic spectrometer. 
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Figure 6: Cross section view of the TWIST experiment indicating the muon beam, time ex¬ 
pansion chambers (TECs), beam scintillator, planar drift and proportional chambers, stopping 
target, and the iron yoke and the magnetic field of the cylindrical superconducting magnet. 
Figure courtesy TWIST collaboration. 


The experiment utilized a 0.7% momentum bite, surface muon beam with a 
near-100% polarization at typical rates of several kHz. To fully characterize and 
periodically monitor the beam properties the beamline incorporated two custom 
time expansion chambers (TECs). Each TEC contained a separate horizontal 
and vertical drift region to measure the horizontal and vertical trajectories of 
individual beam particles. The TECs were constructed with ultra-thin entrance 
and exit windows and filled with low pressure (80 mbar) gas. 

The beam entered the spectrometer through a thin scintillator and a gas 
degrader. A feedback system adjusted the degrader’s helium / CO 2 gas mixture 
to account for environmental changes and maintain centering of the stopping 
distribution in the target foils. 

The magnetic spectrometer permitted the measurement of both the beam 
muons and the decay positron. The setup comprised a symmetric array of pla¬ 
nar drift chambers and planar proportional chambers along the common beam 
and spectrometer axis. The drift chambers enabled the precision determination 
of the momentum-angle coordinates from the helical trajectories of the decay 
positrons. The proportional chambers rendered the \i j e particle discrimi¬ 
nation via the above-threshold time interval of the resulting chamber signals. 
The low-mass chambers were constructed with ultra-thin windows and filled 
with slow-drift gas to optimize position resolution. The support structure was 
constructed from low-thermal expansion ceramics to reduce the effects of tem¬ 
perature fluctuations on spectrometer resolution. 

The experiment utilized two different stopping targets; a 31 jam thickness, 
high-purity silver foil and a 72 jam thickness, high-purity aluminum target. To 
optimize the detector resolution the target foils were employed as the window 
material separating the two central drift chambers of the magnetic spectrometer. 
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A cylindrical superconducting magnet with symmetric upstream and down¬ 
stream holes provided a uniform magnetic field along the beam axis. The field 
uniformity in the tracking region was a few parts-per-thousand and the field sta¬ 
bility was continuously monitored by NMR probes. The field was both mapped 
with a movable Hall probe and calculated using a finite element analysis. The 
track reconstruction used the calculated field distribution anchored by the Hall 
probe measurements. 

The experiment accumulated large datasets of decay positrons from each 
target as well as dedicated measurements of beam properties and other system- 
atics uncertainties. Both experimental data and simulated data were stored in 
a common data format and processed with a common analysis code. 

The analysis involved the reconstruction of the helical particle tracks from 
the raw TDC data and the construction of the positron momentum-angle (p-O) 
distribution from the track parameters. A sophisticated simulation of particle 
interactions, detector geometry, and data readout incorporated: the TEC mea¬ 
surements of beam particle trajectories, the discontinuous nature of chamber 
ionization, a Garfield calculation of the position-dependent drift velocities, an 
OPERA calculation of the magnetic field distribution, and rate-dependent ef¬ 
fects of overlapping tracks. The simulated p-0 coordinates of decay positrons 
included full radiative corrections to order 0 {a‘^) and leading-log radiative cor¬ 
rections to order 0 {a^). 

Importantly the measured and simulated p-^ distributions are functions of 
both the muon decay parameters and the p-0 dependence of the spectrometer 
response function and the track reconstruction efficiency. Therefore the fitting 
of p-0 distributions is used to extract only the differences between the simulated 
values and the measured values of the decay parameters, i.e. Ap, and AS. 
This “relative” method facilitated a blind analysis where the exact values of the 
decay parameters were not known to the collaboration until the completion of 
the data analysis. 

The final results for the Michel parameters extracted from the combined 
analysis of the experimental data and the simulated data were 

p = 0.74977 ± 0.00012 (stat) ± 0.00023 (syst) 

S = 0.75049 ± 0.00021 (stat) ± 0.00073 (syst) 

= 1.00084 ±0.00029 (stat) ± (syst). 

Note that the decay parameter ^ and the muon polarization in pion decay 
are experimentally inseparable. Contrary to p and 5 the uncertainty in 
is dominated by muon depolarization effects including the spectrometer-beam 
alignment and spectrometer fringe fields. 

3.3. Measurement of the positron longitudinal polarization 

A measurement of the positron longitudinal polarization Pl in muon decay 
was recently reported by Prieels et al. [58] . The experimental setup was config¬ 
ured to measure the longitudinal polarization near the positron energy endpoint 
and opposite to the muon spin direction. This kinematical region is uniquely 
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sensitive to the decay parameter that characterizes the angle-energy depen¬ 
dence of Pl. a polarization of 7 ^ 1 for high-energy, backward-angle positrons 
would signal a departure of this decay parameter from the SM value = 1. 

The experiment utilized a high intensity, longitudinally polarized, surface 
muon beam at PSI. The incident muons were stopped in either an aluminum 
target or a sulfur target that were located in a 0.1 T longitudinal holding field. In 
aluminum the initial muon polarization is only weakly depolarized thus yielding 
a large time-averaged muon polarization along the polarimeter axis. In sulfur 
the initial polarization is more strongly depolarized thus yielding a small time- 
averaged muon polarization along the polarimeter axis. 

The detection system consisted of a small acceptance, high precision, cylin¬ 
drical spectrometer followed by a positron polarimeter. The spectrometer was 
orientated to accept only decay positrons that were emitted directly opposite to 
the muon polarization direction. It consisted of three parts: a filter, tracker and 
lens. The filter comprised an initial cylindrical magnet with collimators arranged 
to transmit only high energy positrons. The tracker comprised a superconduct¬ 
ing cylindrical magnet with three planes of position-sensitive silicon detectors 
to measure the trajectories and determine the energies of decay positrons. The 
lens comprised a final cylindrical magnet that focused the decay positrons into 
parallel trajectories at the polarimeter entrance. 

The polarimeter utilized the Bhabha scattering e+e“ ^ e+e“ and in-flight 
annihilation e+e“ ^ 77 of decay positrons on polarized electrons in magnetized 
foils. The cross section for Bhabha scattering and in-flight annihilation is dif¬ 
ferent for the two orientations with the e+e“-spins parallel and the e+e“-spins 
anti-parallel. These spin-dependent cross sections thus induce a signal rate from 
scattering and annihilation that depends on Pl. 

The polarimeter employed (see Fig. a double layer of oppositely magne¬ 
tized foils that were orientated with their magnetization axes at ±45° angles to 
the spectrometer axis. The ±45° configuration yielded an electron polarization 
along the spectrometer axis of roughly 5%. A series of wire chambers and plas¬ 
tic scintillators at the location of the foils was used to discriminate the struck 
foil and determine the scattering / annihilation vertex. A bismuth germanate 
detector array was used to measure the energies and position coordinates of the 
annihilation 7 -rays or the scattered e+e“. Asymmetries in rates of annihila¬ 
tion and scattering were investigated by reversing the foil magnetization and 
inverting the ±45° foil orientation. The measured energies and opening angles 
of 7 -rays or e+e“ pairs—along with the measured energy of the decay positron 
in the magnetic spectrometer—over determined the reaction kinematics and 
allowed for backgrounds suppression. 

The analysis involved constructing the super-ratio 


where is the event ratio from the two foils and the superscript ± denotes 
the two magnetization orientations. The subscript a distinguishes the eight 
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(5) BGO 



Figure 7: Schematic diagram of the double foil positron polarimeter in the Prieels et al. 
experiment m. It shows the two oppositely-magnetized foils, the arrangement of wire cham¬ 
bers and plastic scintillators that discriminate between 7 -ray and e+e“ events and determine 
their vertex coordinates, and the bismuth germanate array. Figure courtesy Prieels et al. 


measurement configurations consisting of (i) ±45-degree foil orientations, (ii) 
aluminum / sulfur targets, and (iii) annihilation / scattering events. The super¬ 
ratios are proportional to the longitudinal polarization Pl of the decay positrons. 
Moreover the use of cancels the effects of differing detection efficiencies for 
scattering / annihilation events from the two foils. 

The key characteristic of differing from unity is the introduction of an 
energy-dependent longitudinal polarization Pl of the backward-emitted decay 
positrons. Consequently, a value of 1 can be identified through a com¬ 

paratively strong energy dependence of for the polarization-preserving Al 
target with a comparatively weak energy dependence of for the polarization- 
destroying S target. No evidence of such an effect was observed in the compar¬ 
ison of the measured super-ratios for the two target materials. 

A combined fit to the super-ratio data for all experimental configurations 
yielded the valu^ 


= 0.981 ± 0.045(stat) ± 0.003(syst). 


®The fitting procedure incorporated energy-independent attenuation factors that empiri¬ 
cally account for globally lower than predicted values for the analyzing powers for the double 
foil arrangement. Although such effects as a smaller than expected magnetization and a 
larger than expected background were investigated, the cause of the attenuator factors was 
not convincingly established. 
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Figure 8: Schematic diagram of the Danneberg et al. experiment m indicating the longitudi¬ 
nally polarized muon beam, beryllium stopping target, transverse precession field, magnetized 
iron-cobalt foil and the bismuth germanate detector array. Figure courtesy Danneberg et al. 


This result represents a ten-fold improvement over earlier experimental work 
and is consistent with the SM prediction = 1.0. 

3.4- Measurement of the positron transverse polarization 

Danneberg et al. [51] have recently reported a measurement of the trans¬ 
verse polarization of the positrons emitted in the /i ^ enn decay. In the limit of 
vanishing positron mass the V-A theory implies neither a T-conserving trans¬ 
verse polarization (z.e. Pti = 0) nor a T-violating transverse polarization (z.e. 
Pt 2 = 0) for the emitted positrons. After accounting for positron mass ef¬ 
fects the average T-conserving transverse polarization in V-A theory is about 
(P^i) = -3 X 10-3. 

The experiment was conducted at PSI. The experimental setup—shown in 
Fig. [S]— incorporated a polarized muon beam, beryllium stopping target and 
positron polarimeter. To optimize the sensitivity to Pt the polarimeter was 
orientated at a 90-degree angle to the muon polarization axis. 

The high-rate muon beam possessed the cyclotron radio frequency time 
structure that consisted of pulses at 20 ns intervals with 4 ns durations. A 
0.37 T transverse magnetic field was used to maintain the ensemble polariza¬ 
tion of multiple muon stops via the matching of the Larmor precession period 
to the 20 ns pulse period. The 4 ns pulse duration of beam particles limited the 
ensemble polarization of stopped muons to roughly 80%. 

The polarimeter employed the detection of y-ray pairs from the in-flight 
e+e- annihilation of the decay positrons on the polarized electrons in a magne¬ 
tized iron-cobalt foil. The resulting angular distribution of gamma-ray pairs is 
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correlated with the relative orientation of the e~^-e~ polarization axes. The time 
and position coordinates of the positron annihilation on the magnetized foil were 
measured by combination of plastic scintillators and planar drift chambers. The 
energy and position coordinates of the coincident 7 -rays were measured using a 
bismuth germanate (BGO) detector array. 

The Larmor precession of stopped muons is important in understanding the 
experimental signature of transverse polarization. This precession would cause 
the corresponding rotation of any non-zero components of the positron trans¬ 
verse polarization. Consequently, the angular distribution of the annihilation 
7 -rays would likewise rotate at the Larmor frequency, with an amplitude and a 
phase that is governed by the transverse polarization {Pti,Pt 2 )- By flipping the 
magnetization of the polarimeter foils or reversing the direction of muon pre¬ 
cession, the associated signal of transverse polarization could be distinguished 
from other systematic effects that might induce periodic signals in the 7 -ray 
time distributions E3 

The data analysis involved constructing time distributions Nij (t) of the 7 -ray 
coincidences in the BGO detector elements (i^j) as a function of the angle ^ be¬ 
tween the 7 -ray pair plane and the horizontal magnetization plane. A non-zero 
transverse polarization would impart on Nij (t) a '^-dependent sinusoidal signal 
with the Larmor frequency and an overall amplitude and phase determined by 
{Pt1,Pt2)‘ 

The experiment recorded about 10^ annihilation 7 -ray pairs from decay 
positrons and permitted measurement of Pti and Pt 2 foi* positron energies 
10-50 MeV. A Monte Carlo simulation was performed to determine the depen¬ 
dence of the transverse polarizations on each decay parameter. A fit was then 
performed of the Monte Carlo polarizations to the measured polarizations to 
determine the “best fit” vaiues for 77 , 77", a'/A and [3'/A. The procedure yieided 
vaiues of 

r} = (71±37±5) X 10“^ 

rj" = (100 ± 52 ± 6) X 10“^ 

for the T-conserving poiarization component Pti and 

olIA = (-3.4 ± 21.3 ± 4.9) X 10“^ 

^jA = (-0.5 ± 7.8 ± 1.8) X 10“^ 

for the T-vioiating poiarization component Pt 2 - Using these best-fit vaiues for 
the decay parameters the authors then derived corresponding vaiues for energy- 
averaged poiarizations of 

(Pn) = (6.3 ± 7.7 ± 3.4) X 10“^ 

{Pt2) = (-3.7 ± 7.7 ± 3.4) X 10“^ 


Specifically, “off-axis” positron annihilation can also induce sinusoidal variations in the 
7-ray time distribution. This off-axis signal was used for polarimeter calibration. 
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The determinations of the transverse polarizations and their decay parame¬ 
ters represent a three-fold improvement over previous experimental work. The 
results are consistent with zero transverse polarization and the SM values for 
the decay parameters, rj = rj" — o' jA — P'IA = {). 

3.5. Global analysis of decay parameters and theoretical implications 

A common representation for the matrix element of the most general, Lorentz 
invariant, derivative-free, current-current interaction is based on definite elec¬ 
tron and muon chiralities 


Moc ^ (15) 

-f=S,V,T 

where T*^, T^, T^ represent the possible (scalar-pseudoscalar, vector-axial, ten¬ 
sor) interactions between a left / right handed muon and a left / right handed 
electron. The parameters are the ten associated complex coupling constants 
of the current-current interaction terms {g^ji and gJ^j^ are identically zero). In 
the case of the V-A weak interaction, the coupling g^^ is unity and all the 
remaining couplings are exactly zero. 

The nine decay parameters from ordinary decay /i+ ^ and one ad¬ 

ditional parameter from radiative decay /i+ ^ e+z/z/ 7 , appear insufficient at 
first glance to unambiguously determine the ten independent complex couplings 
Therefore the global analyses of measured decay parameters involve deter¬ 
mining values for smaller sets of intermediate bilinear combinations of complex 
couplings. The possible limits on non-V-A interactions— i.e. all couplings g]^ 
excepting gY ^—then derived from the values determined for the bilinear 
combinations. 

Sets of intermediate combinations of complex couplings for parameterizing 
the constraints obtained from decay parameters have been introduced by Ki- 
noshita and Sirlin m and Fetscher, Gerber and Johnson m- For example, the 
scheme of Fetscher et al. incorporates the maximum number of positive semidef- 
inite bilinear combinations of complex couplings gf^. Of ten parameters, the 
four bilinears 

Qll = 

Qlr = + \9lr?‘ + 

Qrl = -^\9rl?‘ + \9rl?‘ + 

Qrr = 


^^One coupling is contrained by normalization thus leaving the Fermi constant Gp and 
eighteen additional independent parameters to be experimentally determined. 
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have special significance from their correspondance to the transition probabili¬ 
ties from a muon with handedness /i = L, i? to an electron with handedness 
e = L^R. In the standard model the quantity Qll = ^ and the other bilinears 
are exactly zero. 

In Table we compile the current determinations of the decay parameters 
including the recent results from TWIST, Danneberg et al. and Prieels et al. 
Also listed are the SM values for the decay parameters, i.e. for = 1 and all 
other couplings zero. The world data on decay parameters are consistent with 
the V-A charged current weak interaction to levels of 10“^-10“^ in many cases. 


Table 2: Compilation of the most recent determinations of the nine decay parameters in 
ordinary muon decay and the additional fj decay parameter in radiative muon decay. With 
the exception of the ordinary decay parameter and the radiative decay parameter fj, we 
give the results of the most recent precision measurements. For and fj we give the world 
averages compiled in Ref. | 43| (the average for is dominated by the measurement of Burkard 
et al [60] and the average for r] is dominated by the measurement of Eichenberger et al EU). 


Parameter 

Measured value 

SM value 

Ref. 

p 

0.749 77 ± 0.000 12 ± 0.000 23 

0.75 

m 

5 

0.750 49 ± 0.000 21 ± 0.000 27 

0.75 

[56] 

pfi 

1.00084 ± 0.00029±g:ggJ|| 

1.0 

m 


1.00 ±0.04 

1.0 

m 


0.981 ± 0.045 ± 0.003 

1.0 

EH 

V 

(71 ±37 ±5) X 10-3 

0.0 

m 

rf" 

(100 ±52 ±6) X 10-3 

0.0 

EB 

o! j A 

(-3.4 ±21.3 ±4.9) X 10-3 

0.0 

m 

p'lA 

(-0.5 ±7.8 ±1.8) X 10-3 

0.0 

ED 

V 

0.02 ± 0.08 

0.0 

m 


A global analysis of decay parameters has recently been performed by the 
TWIST collaboration [56]. It followed the procedure of Gagliardi, Tribble and 
Williams [62| but updated their input values for decay parameters with the final 
results of the TWIST experimentThe procedure involved determining the 
joint probability distributions for a hybrid set of nine bilinears of coupling 
from the muon decay parameters via a Monte Carlo procedure (see Ref. [63] for 
details). From these results, they then derived limits on the various couplings 
g]^. Their results for the bilinears Qe/x and the couplings g]^ are reproduced in 
Table [3] 

Of special interest are the upper limits derived on the bilinear sums = 
Qrr Qlr and = Qrr + Qrl A non-zero value of would signal a 
weak interaction contribution from a right-handed muon current and a non-zero 
value of would signal a weak interaction contribution from a right-handed 
electron current. The decay parameter ^ of the positron angular distribution 
is particularily sensitive to and the decay parameter of the positron 


^^This global analysis precedes the recent measurement of by Prieels et al. 
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Table 3: Results for the experimental limits on the bilinear quantities Qe/x and the coupling 
constants of the general current-current interaction derived from the global analysis of 
the muon decay parameters m- In V-A theory the bilinear Qll and coupling are unity 
and all other bilinears and coupling constants are exactly zero. Note the limits on and 

\gY,L\ are derived from inverse muon decay, h'^e~ i^eg~ • 


Qrr < 3.0 X 10-4 

Qlr < 6.3 X 10-4 

Qrl < 0.044 

Qll > 0.955 

|5±I < 0.035 

\91r\ < 0.050 

l5il< 0.420 

\9L\ < 0.550 

\9rr\ < 0.017 

\9lr\ < 0.023 

\9rl\ < 0.105 

\9Il\ > 0.960 

III 

0 

\9lr\ < 0.015 

I^alI < 0.1050 

0 

III 


longitudinal polarization is particularily sensitive to The global analysis 
yielded limits on contributions from right-handed muon interactions of < 
8.2 X 10“^ and right-handed electron interactions of < 0.044. 

Danneberg et al. and Prieels et al. also discuss limits on right-handed cur¬ 
rents involving specific exotic scalar, vector and tensor interactions. For ex¬ 
ample, the decay parameters a'/A and /3'/A derived from the transverse po¬ 
larization Pt 2 are uniquely sensitive to T-violating contributions in purely lep- 
tonic interactions. Danneberg et al obtained a limit on possible T-violating 
scalar interactions between right-handed charged leptons currents of Im(^|^^) = 
(5.2 ±14.0 ±2.4) X 10-^. 

Muon decay is particularly valuable in imposing constraints on various left- 
right symmetric (LRS) extensions of the standard model electroweak interaction 
[64]. Such models introduce a new VpA interaction coupling to right-handed 
currents that partners the known V-A interaction coupling to left-handed cur¬ 
rents in order to restore parity conservation at high energies. In LRS models 
the V-\-A and V-A interactions are mediated by Wr and Wr gauge bosons with 
couplings constants qr and qr- A mass mixing angle ( and CP-violating phase 
uj together determine the relation between the weak eigen-states Wr/r and the 
mass eigen-states W 1/2 


Wr = lFicosC±fF2sinC 
Wr = e^^{-Wi sinC ± W 2 cosC) 

with low-energy parity violation emerging when the kFR-boson mass exceeds the 
kFL-boson mass. In manifest LRS models the two couplings are equal and in 
generalized LRS models the two coupling are distinct. 

A variety of approaches—from direct searches for kFR-boson production 
at pp colliders to setting limits on kFR-boson virtual contributions in K^-K^ 
mixing, ^d-decay and muon decay—have been used in the investigation of the 
various LRS extensions of the standard model (for further details see the re¬ 
views |43l[65]). The results from pp colliders and K^-K^ mixing set impressive 
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Figure 9: Allowed region of W 2 boson mass versus mass mixing angle C in the generalized left- 
right symmetric model. Shown are the limits derived from the recent TWIST measurement 
of Pfj,^ and the earlier results from Beltrami et al. m and Jodidio et al. un- 


lower bounds on ldb?-boson masses while results from unitarity tests that utilize 
0+—nuclear ^d-decay set impressive bounds on Wl-Wr mixing. However, 
unlike these processes the decay of muons is purely leptonic and essentially free 
of any assumptions concerning the CKM matrix elements of the hypothetical 
Wr boson (many analyses assume a Wr boson with standard model-like cou¬ 
plings and CKM matrix elements). Consequently, the TWIST result for [57] 
yields complementary and restrictive limits on the parameter space of the gen¬ 
eralized LRS model. Their limits on QlIsr ^2 versus QlIsr C reproduced 
in Fig. 


4. Charged Lepton Flavor Violating decays 

4 . 1 . Lepton flavor and physics beyond the standard model 

Once a muon, always a muon—or at least the flavor of a muon seems to 
be preserved. The conservation of a “muon” flavor was not expected when 
this particle was discovered, nor does it associate a conserved quantity with a 
fundamental symmetry, as required by Noether’s theorem. But, following non¬ 
observations of its violation over decades of trials in reactions such as the simple 
decay mode /i ^ ey, it emerged as a standard model foundational statement; a 
muon cannot become an electron. 

In a two-neutrino (later, three) world of massless neutrinos, the transfor¬ 
mation of muon flavor to electron flavor (or tau flavor) is strictly forbidden. 
The conservation of separate e— ,/i— , and r-type lepton numbers holds; that 
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Figure 10: a) The standard model allowed decay fi —)■ ey, which proceeds through a loop 
process involving neutrino mixing at the unmeasureably low branching ratio below b) 

A SUSY based diagram depicting smuon-selectron mixing inducing the same process. 


is Al/g = AI/^ = AI/t- = 0 , where the Li are assigned lepton flavor numbers. 
The additive sums are unchanged in any given reaction or decay. The rules 
worked perfectly prior to the period from 1998 - 2001 when the atmospheric, 
solar, and reactor-based neutrino experiments demonstrated that neutrino fla¬ 
vors mix [66j [67l EH] • These discoveries changed the rules. Neutrinos must have 
finite mass, and they can transform within the three generations. In neutral 
processes at least, lepton flavor is not a good quantum numbeip^ 

Violation of lepton flavor in the neutrino sector, in turn, implies that the 
decay /i ^ ey must occur through loop processes. However, using upper lim¬ 
its on the neutrino masses and the measured mixing angles, a branching ratio 
below 10“^^ is deduced for the diagram depicted in Fig. 10 r. Needless to say, 
this is unmeasurably small. The bright side, is that this standard model allowed 
rate—which is representative of the other charged lepton flavor violation (cLFV) 
processes—is so tiny that any non-null measurement must indicate new physics. 
Indeed, calculations based on many popular SM extensions suggest large cLFV 
effects should appear near to current experimental limits. That is the subject of 
the current chapter. The richness and importance of this topic is documented 
in numerous scholarly reviews [saizniiziiizaizH]- Each generally takes on a dif¬ 
ferent emphasis between theoretical expectations, experimental techniques and 
history, or an overall picture. Here, we highlight our view of the key topics be¬ 
ing discussed and the three most sensitive and promising experimental programs 
going forward. 

Comprehensive and high-sensitivity studies of tan lepton cLFV decay modes 
are being carried out now and further plans exist at the upgraded B factory in 
Japan. in. Tagged r decays to muons and electrons, such as r ^ eee, r ^ 
r ^ ey, and r ^ /iy as well as many other combinations are being studied. 
Impressive sensitivity has been achieved with branching ratios typically close to 


current hot topic in neutrino physics is the search for neutrinoless double beta decay, 
If observed, this process violates lepton number conservation and serves to prove that 
neutrinos are their own antiparticles—Majorana fermions. It further supports the leptogenesis 
explanation of why we have a baryon asymmetry in the universe. 
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the 10“^ level and the future promises order-of-magnitude improvements with 
Belle-IL Additionally, the LHCb experiment has already set a competitive limit 
on the BR{t of < 4.6 x 10“^, which was based on the 3fb“^ of data 

acquired in the first LHC run [75]. Presumably this limit will be strengthened 
in the future as the systematics do not dominant the uncertainty. 

However, it is the three low-energy muon reactions 

1. /i+ ^ e +7 {BR :< 5.7 x Ref. |8] 

2. /i+ ^ e+e“e+ {BR :< 1.0 x 10“^^) Ref. |9] 

3. fi- ^ e-N {BR :< 7 x Ref. [lO] 


that have the more impressive limits and greater sensitivity to new physics. 
Processes 1 and 2 represent ultra-rare decay modes. Process 3 involves the 
formation of a muonic atom with a nucleus—the limit quoted here is from the 
/i“ An atom—followed by coherent conversion of the muon to an electron, which 
is ejected with an energy close to the muon rest mass. The history of these 
measurement sensitivities is compiled in Fig. 11, which includes projections for 
experiments now being upgraded or constructed. 

Note that other processes also exist. For example, the spontaneous con- 


muonium atom to anti-muonium, fi e+, has been 


version of the M = /i+e“ 
searched for and a limit of Pmm ^ x 10“^^ (90% C.L.) in a 0.1 T magnetic 
field has been set [76]. Here, electron and muon number are violated in the 
same exotic process. The physics reach is not as competitive if directly com¬ 
pared to current cLFV efforts and there are no current plans to improve it on 
the horizon m- But, it should be noted that this is a rather complementary 
process that might be induced by different physics; for example, see m 

A relative comparison remark is in order here. Generically, for a loop-induced 
process, similar to those that might cause the deviation from the SM for the 
muon anomaly, the expected rates for processes 1:2:3 scale as 389 : 2.3 : 1; 
(here, fiAl is assumed for 3) [70]. The ratio implies that the present /i+ ^ e+y 
limits greatly exceed current limits from the 3e decay or the fi ^ e conversion. 
However, the ambitious goals of modern fi ^ e conversion experiments aim for 
4 orders of magnitude improvements, which will close the gap. 

The comparison using loop-like scaling is not the whole story, since the 
sensitivities vary depending on what kind of new physics leads to cLFV. Ulti¬ 
mately, one would like to have measured deviations in all three channels to sort 
out the underlying mechanism. Although there are many compelling standard 
model extensions, such as the SUSY process involving smuon-selectron mixing 
as depicted in Fig. [Tq]3 , we will follow the more generic approach developed by 
de Gouvea HU. In this model-independent analysis, the cLFV processes can be 
compared for their respective new-physics sensitivities in terms of an effective 
energy scale versus a parameter that slides from loop-like exchanges to contact 
interactions, de Gouvea considers an effective Lagrangian to describe processes 
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Figure 11: The solid shapes indicate the achieved 90% sensitivity limits for ^ e +7 
(triangle), —)■ e+e“e+ (circle), and the coherent conversion process ii~N e~N (square). 

Expected improvements from approved experimental programs have light grey interiors and 
the same shapes. Figure updated from EH, courtesy R. Bernstein. 


(1) and (3) of the form 

^cLFV = 7- , 

[K + 1 )A^ 

+ 7^1)A2 PL7fieL{uLj^UL + 

Here L and R are fermion field chiralities and is the photon field. If the 
dimensionless parameter k 1 the process is dominated by a magnetic moment 
type operator (as above). In contrast, if k ^ the four-fermion operators (2nd 
line) dominate, representing a point-like contact interaction. The parameter 
A is an effective energy scale that extends for current experimental sensitivity 
goals to lOOO’s of TeV—for optimized coupling—well in excess of any collider 
reach. This comparison is plotted in the left panel of Fig. [^for both current 
experimental limits and the future projects described below. The sensitivity 
from /i+ ^ e +7 searches falls off quickly for contact-like interactions, unlike 
fi — e conversion, which rises. 

A similar Lagrangian (the second term is different) can be used to describe 
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reactions ( 1 ) and ( 2 ): 


^cLFV = (k+ 

where n and A correspond to similar definitions. The right panel of Fig, plots 
these sensitivities. The fi eee process will be several orders of magnitude less 
sensitive compared to ^ e +7 for loop-like processes; however, for contact- 
type interactions there is a part of the new-physics search space that favors the 
3e decay mode. 

An important word of caution is in order here so as not to over-interpret 
the vertical axis in Fig. The new physics “mass scale” Aclfv is not to be 
interpreted strictly as mass. It is product of mass and flavor-violation mixing, 
where the latter is similar to the mixing angles in a CKM or PMNS matrix for 
quarks and neutrinos. The mixing can be large—neutrino like—or very small— 
quark like, or it can be anything at all. With a cLFV experiment there is but 
one-measurement and two parameters, which emphasizes the need for multiple 
measurements (channels) and comparisons to other precision measurements. 
For example, consider the situation where the new physics impacts both g — 2 
(flavor conserving) and cLFV. Then the new physics scale can be related by 
©e^A^Lpy = A^_ 2 , where ©g^ represents the level of flavor-mixing in the new- 
physics process; ^ 1 is a flavor indifferent, maximally mixed scenario, and 
©e^ ^ 1 implies small mixing m- 


4 . 2 . Experimental ehallenge of muon eLFV 

The goals for modern charged lepton flavor violation experiments all aim 
at single event branching ratios of a few times 10“^^ to 10“^^ range. To be 
sensitive to such stunningly rare processes requires an event signature that is 
unique with respect to backgrounds, and exceptional design and execution by 
the experimental team. Because all processes involve muons at rest, the energy 
scale of the emitted particles is low—below for the decay studies and equal 

to ^ for the conversion process. From a detector design perspective, this is a 
challenge, requiring ultra-thin tracking detectors and especially high-resolution 
calorimetry. The sub-dominant, but ordinary, muon decay modes are problem¬ 
atic. Radiative muon decay, /i+ ^ occurs with a branching ratio of 

1.4(4)% and radiative decay with internal conversion, ^ e+e“e+z/gZ/^, oc¬ 
curs at the (3.4 ± 0.4) x 10“^ level. For the corners of phase space where the 
neutrinos escape with very little energy, these decay modes can mimic the cLFV 
processes ( 1 ) and ( 2 ), at a many orders of magnitude greater rate. 

To achieve sensitivities below 10“^^ a considerable number of stopped muons 
is required (well above 10^^ after factoring in efficiencies). In context, the MuLan 
lifetime experiment described in Sec. 2 . 2.1 measured to 1 ppm using 2 x 
10^^ decays, and utilizing a high-intensity beamline at PSI. The cLFV program 
requires, in contrast, muon samples larger by factors of 10^ — 10^. The highest- 
intensity 7 rF ^5 beamline at PSI can deliver a dc rate of 10 ^ muons/sec, which 
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Figure 12 : Comparison of the new physics mass scale reach for the cLFV experiments 

using current achieved sensitivities and goals of planned programs. In both panels, k ^ 1 
corresponds to loop-like processes, transitioning to 4-fermion contact interactions for large 
values of k. Left panel: —)■ e +7 versus iJ.~ N e~ N coherent conversion experiments. 

Right panel: —)■ e +7 versus e+e“e+. Figure updated from ED, courtesy A. de 

Gouvea. 


is 10 times higher than MuLan used. The dc nature of the beam structure 
is ideal for the decay programs that observe one event at a time. However, 
at high rates, multiple events are in the detector at once, requiring excellent 
coincident timing and position resolution to distinguish individual decays. To 
reach the ultimate goal of the 3e experiment of 10“^^ will require a new ultra- 
high intensity beamline having 10 x greater flux than presently exists. For the 
next-generation /i ^ e conversion experiments—the goals there approach the 
10“^^ level—an entirely different concept is required, which we describe below. 


4 . 2 . 1 . Measuring /i+ ^ 
The decay mode 


e+7 


e +7 is being pursued by the MEG collaboration at 


PSI. Their results to date set the limit < 5.7 x 10“^^; 90% C.L. [8]. The 

decay mode features the back-to-back emission of a positron and a gamma ray, 
each having energies equal to m^/2 ^ 53MeV. Positive muons at a dc rate of 
up to 30 MHz are stopped in a thin polyethylene disk centered in a cylindrical 
geometry. A tracking drift chamber system and a high-resolution liquid xenon 


(LXe) calorimeter are key components, see Fig. 13 


A unique feature of the experiment is the C Oust ant Bending RAdius (CO¬ 
BRA) magnet. Its gradient field shape is designed to maintain a near-constant 
bending radius of decay positrons, independent of their initial pitch angle, and 
further, to sweep away the decay trajectories to the upstream or downstream 
side of the stopping target. The unique calorimeter uses 862 PMTs, submerged 
in LXe and mounted on all surfaces to sum the light and to therefore provide en¬ 
ergy, time-of-arrival, and location of the interaction of the gamma in the xenon 
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Figure 13: The side and front schematic of the MEG detector with a hypothetical /r —)■ ey 
event superimposed. Figure courtesy MEG collaboration. 


bath. The ultra-low-mass drift chamber sits in a half circle below the target. 
A scintillator array is used to precisely establish the timing needed. At signal 
energies = Eq+ = m^l2^ the tracker resolution is 1.5% and the calorimeter 
resolution is 4.5%. The overall acceptance for decays of interest is ^ 18%. [78] 

Two types of backgrounds enter, those from radiative decay, and those from 
pileup events. The key to distinguish these from signal are the performance 
parameters of the detectors. We reproduce Fig. 14 from [ 8 ] because it perfectly 
illustrates the case by displaying the allowed signal event regions, which are 
dictated by detector resolutions. If the resolutions can be improved, the signal 
region can be reduced, increasing the overall sensitivity. The right panel shows 
the cos&ej vs. At plane, where ©e^ is the difference between the emitted angles 
of the electron and gamma, and te^ is the difference in their timing. The left 
panel shows the e+ — 7 energy plane with the enclosed ellipses defining the signal 
region. Improvements in energy-, angle-, and timing-resolution of the involved 
detectors allows one to shrink the good-event windows proportionally. 

The MEG II approved upgrade program m is designed to handle muon 
stopping rates on a thinner target at rates up to 70 MHz; the final statistics will 
require three years of running, beginning in 2016 . Resolutions will be improved 
by about a factor of 2 on all detectors. New UV-sensitive SiPMs will replace the 
PMTs and provide more uniform light collection. A single-volume drift chamber 
system will replace the current vane structures. The positron scintillator coun¬ 
ters provide 30 ps timing resolution. The overall goal for the branching ratio is 
< 4 X 10 -^^. 


4 . 2 . 2 . Measuring e+e“e+ 

The unique signature of the e+e“e+ reaction is the three coincident 

e^s emanating from a common vertex and carrying, in sum, a total energy cor- 
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Figure 14: The event search regions from the MEG publication [ 8 ], which illustrates the 
importance of a comprehensive approach to improving detector resolutions and timing in 
order to narrow the allowed good event region. The left panel shows the — Ee plane, where 
the allowed region is centered at m ;^/2 in both projections. The right panel is the At— cos ©67 
plane where the region of interest picks out simultaneous, back-to-back decay products. The 
blue enclosed regions are signal probability distribution function contours of 1 — , 1.64—, and 
2 ( 7 . Figures courtesy MEG collaboration. Ref. [H] 

responding to the full muon mass, and a net momentum of zero. This measure¬ 
ment is a particular challenge, as the three tracks must be accurately measured 
with very high resolution for each event, and the triple coincidence implies that 
the tracker efficiency be very high. Candidate three events having a to¬ 
tal energy near to the muon mass must be considered. Background processes 
such as the internal-conversion decay, generally yield lower event energy sums 
and distort the momentum balance such that the event should be outside the 
search window. Multiple scattering and energy loss in the tracking chambers 
could distort the interpreted energy if the detector design does not factor that 
in properly. Of the three processes discussed, this one has been idle for a long 
time—see [9] with BR^e < 10“^^—, but the new Mu3e collaboration at PSI 
is seeking to study it in stages. They are approved with a goal of 10“^^ using 
existing beamlines. A final phase is imagined but requires a muon beam with a 
rate of 10^ Hz, which is a non-trivial development and investment [80] . 

The key to the ultimate precision is an extremely high muon stopping rate 
distributed over a double cone thin aluminum target. The decay trajectories are 
then imaged precisely using state-of-the-art thin silicon pixel detectors assem¬ 
bled into cylindrical geometries. Timing scintillators between tracking layers 
mark the event times. To achieve the ultimate sensitivity of 10“^^, the resolu¬ 
tion on the reconstructed muon mass must be roughly 0.5 MeV (the mass of a 
single electron!), just to suppress the internal conversion look-a-like decay mode 
/i+ ^ This further assumes excellent sub-ns coincident timing of 

the tracks, and does not factor in accidental backgrounds. 

A schematic diagram of the approved Mu3e experiment is shown in Fig. \TE\ 
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Figure 15: Side and end views of the proposed Mu3e detector with a single simulated 

—)■ e+e“e+ event displayed. The conical target distributes the high rate of muon stops 
in different locations forming distinct vertices in the reconstruction of tracks. The ultra- 
thin pixelated silicon tracker extends up- and down-stream. The scintillating fiber hodoscope 
provides precision timing. The experiment will be assembled and run in phases; the approved 
Phase IB configuration is shown. Figure courtesy Mu3e collaboration. 


In the first phase, only the central pixel detector surrounding the target will be 
used and the beam rate in the nEb line will be limited to 10^/i+/s. The physics 
goal is a BR sensitivity of roughly 10“^^. The next step increases the tracking 
volume, adds scintillating-fiber trackers and other components, and increases 
the muon stopping rate to 10^ Hz. The aim is a sensitivity of 10“^^. 


4 . 2 . 3 . Measuring ji N ^ e N coherent conversion 

The quantity of interest in a /i ^ e experiment is the ratio 


, ^ r[/i- + A{Z, N)^e- ^ A{Z, N)] 

~ r[/i- + AiZ, N)^i^^ + A{Z - 1, N)] 


(16) 


where the rate of ordinary capture in the denominator is already quite well 
known. For example, for typical stopping targets such as aluminum, which we 
will continue to assume below, one can observe that the muon lifetime of a /i-Al 
atom is reduced to ~ 864 ns. From 

~ Ttot ~ Tdecay T Tcapture 


and knowing the free lifetime, one obtains the capture rate T^ai ~ 7 x 10 ^ s“^. 
The numerator in Eq. on the other hand, takes on a new level of challenge. 

While the history of /i ^ e experiments is well covered in Ref. [72], here we 
concentrate on the two major effort ^^t hat ambitiously aim at single-event sensi- 


^"^We note that the less ambitious DeeMe m experiment at the J-PARC Material Life 
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tivities approaching ^ 3x 10“^^, a 4 order of magnitude improvement compared 
to present limits. The two projects—Mu2e [82] at Fermilab and COMET [83] at 
J-PARC—share many similar conceptual design features, while individual tech¬ 
nical solutions differ. We discuss the generic method here using some examples 
for clarity. 

A principle requirement is a sample of ^ 6 x 10^^ muonic atoms. No ex¬ 
isting secondary decay beamline can provide such a flux to accumulate that 
many events in a few years of running. Furthermore, the muons should be 
delivered with a beam-on / beam-off time structure commensurate with the 
muon lifetime in the target material. Ideally, bursts of muons should arrive and 
stop in a suitable nuclear target during an accumulation period lasting some 
lOO’s of ns at most, followed by a beam-off and background-free quiet measur¬ 
ing period. During that time, the detector is sensitive to the unique /i ^ e 
signature event having a mono-energetic electron emitted with Ee = 105 MeV; 
the value represents the muon mass, less the atomic binding energy. For maxi¬ 
mally efficient data taking, the combined accumulation-measuring cycle should 
be 0(2r(^y^p) ~ 1.7/is and repeat as efficiently and continuously as possible, 
subject to the main accelerator macro cycle timing. 

To provide some numbers as an example, the Mu2e statistical budget as¬ 
sumes an average data-taking rate of ^ 10^^ formed muonic atoms per second 
and 6 x 10^ s of running time. Using the reduced fraction of the Fermilab 1.33 s 
accelerator cycle available, they must collect nearly 60,000 atoms in a 1.7/is 
cycle. How is this possible? 

Both Mu2e and COMET follow a recipe credited to Lobashev and proposed 
for MELC in 1989 [84]. A pulsed proton beam is slammed through a target that 
is enveloped in a superconducting “production” solenoid. The produced pions 
are contained by the magnetic field, as are their decay muons below a momentum 
of ^ 40 MeV/c. The solenoid field has a strong gradient to direct the spiraling 
secondaries out along the upstream direction into a curved “transport” solenoid. 
This either S— (Mu2e) or U— (COMET) shaped device provides momentum and 
sign selection of the captured muons and rejects by line-of-sight the transport 
of neutral particles from the production target. It is long enough to allow most 
of the pions to decay before emerging at the exit. The final stage is yet another 
solenoid—the “detector” region—where the muons will first stop in thin A1 
blades. In the same volume, and downstream further, a specialized spectrometer 
system is positioned to intercept the highest energy decay electrons (well above 
the Michel endpoint). Figure 16 illustrates the geometry and scale of the Mu2e 
interpretation of this concept. As noted, the transport section for COMET is 
U shaped,and the detector section is also U shaped, the bend between target 
and tracker. 

Assuming now that one has a sample of formed fiAl atoms, with the muon 


Science Facility aims to turn on sooner than these efforts, with the single event sensitivity 
of 2 X 10“If achieved, that level will be more than an order of magnitude improvement 
compared to current limits. 
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Figure 16: The Mu2e Experimental concept. The pulsed proton beam enters the production 
solenoid and passes through a tungsten rod. Pions produced there, decay to muons which are 
directed via magnetic gradients into the S'—shaped transport solenoid, arriving in the detector 
solenoid region. The momentum-selected negative muons stop in an A1 target. The signature 
decay 105 MeV electron is measured with a tracking chamber followed by an electromagnetic 
calorimeter. Figure courtesy Mu2e collaboration. 


having cascaded to the IS ground state, the next critical issue is associating 
a high-energy emerging electron with /i ^ e conversion. One might naively 
assume this is quite far away from the Michel endpoint which 

it is; however, that only applies for an unbound muon. The muon in a /iAl 
atom can “decay in orbit” (DIO) about 39% of the time. The corresponding 
kinematics of the in-motion system allows—in principle—the emitted electron 
to carry away the entire muon mass in energy, with the nucleus participating 
to conserve momentum. This is a well-known effect and the probability of 
an electron close to the 105 MeV signal region is indeed quite rare (or else, 
conversion experiments would have long ago ceased). With the approaching 
efforts of Mu2e and COMET, the exactness of the DIO spectrum had to be 
revisited to determine the real sensitivity of these major new efforts. In a 
recent calculation, Czarnecki, Tormo and Marciano included proper nuclear 
recoil effects and generated the expected spectra of electron rate vs. energy 
shown in Fig. The left figure provides the full range of energies; the right 
panel is a blowup of the critical region from 100 to 105 MeV. The scale is in 
1 MeV bins of the relative sensitivity units projected by the experimentalist. It 
is clear that very high resolution is thus required above 100 MeV in order to 
keep the DIO fraction below the single event sensitivity. At the time of this 
writing, both Mu2e and COMET are in construction phases. 


5. Muon Dipole Moments 

5.1. Terminology for muon magnetic and electric dipole moments 

The muon’s magnetic dipole moment (MDM) is related to its intrinsic 
spin s with a proportionality that includes the ^-factor, which is embedded in 
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Figure 17: Left: Decay in orbit electron energy spectrum assuming an aluminum target 
for the /rAl atom. Right: Blowup of the critical signal region from 100 to 105 MeV where 
the black squares include recoil effects; red triangles do not. Figures courtesy A. Czarnecki, 
Ref. [85] 


the relation 


(&;)*■ 


( 18 ) 


The provision for a possible nonzero electric dipole moment (EDM) with a 
magnitude parameterized by 77 , is given by 


(^) 


(19) 


Notice that any EDM must be aligned along the angular momentum axis, the 
only vector in the system. We first discuss the important topic of the MDM, 
where the measured value is in disagreement with the standard model (SM) 
prediction. The measured EDM is compatible with zero, but it is not yet de¬ 
termined to high enough precision to challenge SM completeness. The EDM 


prospects are discussed in Sec. 5.4 


5.2. Magnetic dipole moment 

The Dirac equation predicts = 2 for the structureless, spin- 1/2 muon. 
Radiative corrections from electromagnetic (QED), weak and hadronic loops 
give rise to a so-called anomalous magnetic moment a^; that is, g^ = 2(1 + a^), 
or more commonly stated: = {g — 2)/2. The anomaly is a small correc¬ 

tion, with ~ 1/850. Eigure 18 illustrates example Eeynman diagrams that 
must be evaluated to arrive at the SM value for a^. In practice, more than 
10,000 of these topologies have been calculated. These include QED through 
5 loops (10^^ order), weak exchanges through 3 loops, leading-order hadronic 
vacuum polarization (Had-LO)—which is determined from experimental data— 
and higher-order hadronic light-by-light scattering (HLbL), which is difficult to 
evaluate owing to the non-perturbative nature of QCD at low energies. An 
excellent and detailed summary can be found in the textbook and review by 
Jegerlehner [86l|87|. Eigure [^displays the magnitudes and uncertainties of the 
SM contributions along with the current summary. 
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Figure 18: Example diagrams that contribute to the muon anomalous magnetic moment, a) 
Leading-order QED “Schwinger” term; b) Electroweak Z exchange diagram; c) Lowest-order 
hadronic vacuum polarization; d) Hadronic light-by-light scattering. 


The experimental history of {g — 2) measurements is quite rich and well 
reviewed; see, for example Ref. [88]. A series of storage ring experiments was 
conducted at CERN in the 1960’s and 70’s and the highest precision effort was 
completed at Brookhaven National Laboratory in 2001. As this review will 
describe, next-generation efforts are being prepared niEn]. The vertical black 
bars in Figure [T^ indicate the achieved uncertainties of the three most recent 
experiments (solid lines) and the precision goal (dashed line) of the new Fermilab 
experiment. 

Both experiment and theory are now known to similar sub-ppm uncertainty 
and the comparison provides a sensitive test of the completeness of the standard 
model. If the SM accounting is accurate, and the experimental result is correct, 
the present comparison already begins to suggest the existence of some new 
physics process that affects a^. 

Currently accepted SM predictions for a^(SM) are 

a^(SMa) = 1165 918 02(49) X 10““ (0.42ppm) and 
a^(SMb) = 1165 918 28(50) X 10““ (0.43ppni). 

The subscripts a m) and b ([92|) represent slightly different evaluations of 
the leading-order hadronic vacuum polarization contributions. In both cases a 
theoretically sound dispersion relationship is used, which evaluates from an 
appropriately weighted integral of e+e“ ^ hadrons absolute cross section data, 
summed over all energies. The uncertainty here is largely experimental. The 
HLbL contribution, on the other hand, is model based and its quoted theoretical 
uncertainty is only estimated at present. 

The experimental is based solely on the BNL E821 measurement |6|, 

a^(Exp) = 1165 92091(63) x 10“^^ (0.54ppm). 

The difference between experiment and theory—the “(^ — 2) test”—is 

Aa^(Exp — SMa) = 289(80) X 10“^^ (3.6 a) or 
Aa^(Exp - SMb) = 263(80) x 10“^^ (3.3a). 
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Figure 19: The horizontal bars display the magnitudes of the standard model contributions 
(blue) and their present uncertainties (red). Terms include QED through 10^^ order, leading- 
and higher-order weak, leading- and higher-order hadronic vacuum polarization, and hadronic 
light-by-light. The 2014 summary is given along with its uncertainty of 49 X 10“^^. The 
solid vertical black lines are placed at the achieved final precisions of the CERN II (27-ppm), 
CERN III (7.3 ppm), and BNLE821 experiments (0.54 ppm). The dashed black line illustrates 
the precision goal of FNALE989 (0.14 ppm) 
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In both cases, the result is quite provocative. The magnitude of the difference is 
large—it exceeds the electroweak contributions—and the statistical significance 
is large. This persistent discrepancy has led to many speculations of new physics 
scenarios and to challenges to the accounting procedures, data selection, and 
hadronic models used to determine the SM expectation. The resolution—new 
physics or some type of theory or experimental error—has led to the launching 
of two major new experimental thrusts, with one of them that aims to reduce 
(5a^(exp) fourfold to an absolute uncertainty of 16x 10“^^. In parallel, a vigorous 
theory campaign is taking place [93l [94] . 

As Fig. [^illustrates, the QED and weak uncertainties are already well below 
any experimental reach. Improvements by a factor of 2 in the dominant Had-LO 
uncertainty are expected from new data sets at VEPP-2000 (Novosibirsk) and 
BESIII (Beijing), along with continued analyses of BaBar and Belle existing 
data sets [93] . 

While the magnitude of the complete HLbL contribution is ^ 110 x 10“^^, 
the quoted uncertainties can be as large as 35%. A recent review of the status 
of the HLbL contributions is given in [95]. Initiatives going forward include 
those based on lattice QCD [96], and new 7 * physics measurement programs 
at BESIII and KLOE (Erascati) that aim to build a data-driven approach to 
leading HLbL terms [93] . 

5.2.1. New physics possibilities from {g — 2 ) 

What is nature trying to tell us if the current discrepancy Aa^ remains 
as large as it is and the significance eventually exceeds 5 a? While numerous 
explanations exist in the literature, it is instructive to take a more generic ap¬ 
proach here as first outlined by Czarnecki and Marciano m, and elaborated 
on by Stockinger [98]. The magnetic moment is a flavor- and CP-conserving, 
chirality-flipping, and loop-induced quantity. Any new physics (N.P.) contribu¬ 
tions will typically contribute to as 


Sa^{N.P.) = 0[C{N.P.)] X 



( 20 ) 


where M is a new physics mass scale and C{N.P.) is the model’s coupling 
strength, which is then also common to the same new physics contributions to 
the muon mass. That means 


C{N.P.) = 


Sm^iN.P.) 


rriu 


( 21 ) 


Different predictions for a^{N.P.) are illustrated in Eig. 


20 for various con- 
limits, the latter 


pling strengths, C{N.P.) versus the present and future 
being the combined uncertainty from SM theory and experimental sensitiv¬ 
ity. Eor radiative muon mass generation, C{N.P.) = 0(1), which implies 
probes the multi-TeV scale. Eor models with typical weak interaction coupling, 
C{N.P.) = O^ajAn)^ the implied mass scale is very light, arguably ruled out 
by direct measurements. In contrast, models with enhanced coupling such as 
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supersymmetry, unparticles and extra dimensions are represented by the central 
band, where the overlap with corresponds to the TeV-scale physics regime. 
Focussing on SUSY, the expected contribution to has the following behavior: 

of S’- ^ 130 X I0-“ tmUsimip). (22) 


Here, tan /3 is the ratio of vacuum expectation values for the two Higgs doublets 
and sign(/i) is the sign of the higgsino mass parameter, taken to be +1 from 
current {g — 2) constraints. A wide, yet natural, range for tan /3 of 5 —50 provides 
the width of the central band. It is clear from this figure that the LHC reach 
overlaps well with the simplest SUSY expectations, which have been highly 
motivated given the compatibility with the current {g — 2) result [99]. However, 
the lack of any signal in the 7-8 TeV data taking at the LHC has pushed the 
mass scale near to and even above 1 TeV, which has also spawned many variants 
on the simplest SUSY models that must remain compatible with a wide suite 
of experimental limits and the {g — 2) “signal;” see reviews: |1QQ[ llQlj . 

Hypothetical dark photons are very weakly interacting and very light parti¬ 
cles that could produce a large enough contribution to to explain the discrep¬ 
ancy between experiment and theory |lQ2j . In Fig. 20, a dark photon mechanism 
would correspond to a new band to the left of the red band that crosses the 
Aa^ region in the 10 - 100 MeV mass range. The coupling strength would be 
appropriately tuned and small. While initially suggestive of a neat explana¬ 
tion for {g — 2) recent experiments |103j have all but ruled out the simplest 
versions of the theory and its implied parameter space; however, more complex 
scenarios—ones that imply a dark Z' for instance—remain viable m- 


5.3. Experimental eonsiderations 

Given the high impact of the {g — 2) experiments, numerous general and 
detailed reviews have been written [1061110711108j , along with a comprehensive 
publication by the E821 collaboration [6]. Here we will describe the essential 
features that enable such a precision measurement. The persistent > 3cr dis¬ 
crepancy between experimental and theoretical over the past 10 years has led 
to the development of two new experiments that are being designed to resolve 
the situation. FNAL E989 [89] and J-PARC E34 [90] are following different ap¬ 
proaches. In both cases, intense bunches of polarized muons are injected into a 
uniform-field magnet, where the muon spin precession frequency is proportional 
to a^. The decay positrons are measured as a function of time, the rate of which 
features a sinusoidal modulation imprinted on the exponential decay. Beyond 
this, the difference in implementation is quite significant. We first outline com¬ 
mon considerations for a {g — 2) experiment and follow with brief descriptions 
of the two new experimental campaigns. 

The spin of a muon at rest will precess in a magnetic field B at the Larmor 
frequency, ujl = —g(lBl2m^. A precision measurement oi ujl together with an 
equally precise determination of B gives g. Since {g — 2) is known to sub-ppm 
already, a direct measurement of g would need to be at the sub-ppb level to be 
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Figure 20: Left: Generic classification of mass scales vs. contributions from new physics 
sources. Green line: radiative muon mass generation; Red line: Z'^W', universal extra 
dimensions, or Littlest Higgs models with typical weak-interaction scale coupling; Purple 
band: unparticles, various extra dimension models, or SUSY models where the coupling is 
enhanced. The width illustrates a tan /3 range of 5—50 for SUSY models. The yellow horizontal 
band corresponds to the current difference between experiment and theory and the blue band 
is an improvement with a combined theory and experimental error of 34 x 10“^^. Figure 
courtesy D. Stockinger; see also [105]. Right: Neutralino-smuon (a) and Ghargino-sneutrino 
(b) one-loop SUSY contributions to a/^ m 


competitive. However, the muon mass is “only” known to 34ppb, so that 
quantity too would need a major improvement. 

In contrast, a measurement using in-flight muons in a magnet can directly 
determine {g — 2). It is based on the difference between the cyclotron and spin- 
precession frequencies for a polarized ensemble of muons that circulates in the 
horizontal plane of a uniform vertical magnetic field. The cyclotron frequency 
when H = 0 is 


and the spin turns at frequency 


2 m ym 

owing to the torque on the magnetic moment and including the Thomas pre¬ 
cession effect for the rotating reference frame |lQ9j . 

If g were exactly equal to 2 , Ug = however, for ^ 7 ^ 2 , 


CJc = -- 


my 


(23) 


fflB 




qB 






fl - 

2 ) m m ’ 


(25) 


where we have defined uja as the anomalous precession frequency. It is this 
quantity that must be measured to determine {g — 2). 
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The recent storage ring experiments used electric quadrupoles to provide 
vertical containment—effectively creating a large Penning trap. The motional 
magnetic field seen by a relativistic muon in an electric field E will contribute 
an important term to the spin precession rate. Additionally, a nonzero muon 
EDM will also require modification to Eqn. The full expression is then 


^net 

m 



1 \ /? X E 

7 ^ — 1 / c 


n 

2 



(26) 


where iUnet = E Oedm- 

The proposed J-PARC E34 experiment will not use electric focusing, which 
simplifies Eqn. and separates the contributions of uja and uoedm- In contrast, 
the modern storage ring experiments operate with a relativistic gamma of 29.3, 
(P^ = 3.094 GeV/c), which makes the first expression in parentheses vanish. In 
practice, a combined correction to uJa in BNLE821 owing to the electric-field 
and related pitch correction was +0.77 ± 0.06 ppm 

Parity violation in the muon decay chain /r e provides the neces¬ 

sary polarimetry that is required to access the average muon spin direction vs. 
time; that is, the link to uOa- The CM correlation between the emitted angle 
and energy of the decay electron with respect to the muon spin direction is illus¬ 
trated in Eig. To a good approximation, the electron energy in the boosted 
laboratory frame is related to the energy and angular distribution in the CM as 

Ee,lab ~ + COS^cm)- (27) 

This is the key relation. The highest energy electrons are preferentially emitted 
when the muon spin is opposite to its momentum (alternatively, along its mo¬ 
mentum for e+ from /i+ decay). Every decay electron has a momentum smaller 
than its parent muon and consequently curls to the inside of the storage radius 
where a detector is positioned to intercept it and measure its energy and arrival 
time. 

The accumulated number of electrons having an energy greater than Eth 
forms a distribution vs. time that has the structure 

N{t; Eth) = [1 + ^ cos{uJa,t + </>)] (28) 

for a 100% polarized beam. The normalization Nq and asymmetry A depend 
on Eth. The ensemble-averaged spin direction at t = 0 is represented by 0, 
which can have a subtle energy dependence because the time of decay measured 
at the detector might have an energy-dependent time-of-flight component from 
the time of the muon decay. A representative data set from BNLE821 is shown 
in Eig. [ 21 ^. The time-dilated lifetime of ^ 64.4 fis is evident, upon which 
is the modulation from The actual N{t;Eth) distribution can be more 


^^The finite momentum spread, ~ 0.15%, means the term in parenthesis does not 

vanish completely. The vertical betatron oscillations cause B • P^ to not always equal zero. 
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complicated compared to Eqn. [^because of coherent betatron oscillations that 
give rise to further modulations of A and The statistical uncertainty 
on uja has been described in detail in Ref. mni. Different fitting methods and 
weighting schemes can produce reduced statistical uncertainties from the same 
data; however, sensitivity to leading systematic errors typically increases when 
one is using the more aggressive analysis procedures. The simple threshold 
method that is robust and tested has a relative uncertainty on SuJal^a that 
behaves as 

\l NA^P)^^ 

where we include the ensemble averaged polarization (P) for completeness. To 
minimize dujal^a^ h is advantageous to: 1) Use a high magnetic field {uJa oc P); 
2) Run at high energy (increases yr^); 3) Employ a highly polarized muon 
source; 4) Optimize the figure-of-merit (EOM), N. The latter occurs when 
EthlEmax - 0 . 6 . 

To obtain from the measurement of uJa requires an equally precise mea¬ 
surement of the magnetic field. This was accomplished in BNLE821 using a 
suite of pulsed NMR probes. They were used to establish the absolute field 
magnitude, control the time stability of the field via feedback to the magnet 
power supply, and to periodically map the field in situ by use of a multi-probe 
NMR trolley that could traverse the circumference of the storage ring orbit with¬ 
out breaking the vacuum. The water or petroleum jelly filled probes provided 
a field value in terms of the free proton precession frequency ujp. Minimization 
of the systematic uncertainty in the field integral is intimately related to the 
intrinsic field uniformity. Minimizing multipoles having higher order than the 
dipole depends on the care and precision of the shimming tasks carried out 
prior to physics data taking. The final azimuthally averaged magnetic field is 
represented by a contour map as shown in Eig.|2l]3. 

With uja and ujp measured, the muon anomaly was obtained from 


Wa/ujp 


R 


I 


X-R 


(30) 


In the ratio R = Wa/cJp, cOp is the free proton precession frequency in the av- 
erage magnetic field experienced by the muons. The muon-to-proton magnetic 
moment ratio A = 3.183 345107 (84) is determinecj^ from muonium hyperfine 
level structure measurements [m] together with QED, see Sec.|^and Ref.g]. 


5.3.1. Fermilah muon {g — 2) experiment 

The precision goal of E989 is 5a^ = 16 x 10“^^ (140ppb), a four-fold im¬ 
provement compared to BNLE821. This error includes a 100 ppb statistical 


^®True for the storage ring experiments that use quadrupole focussing; absent for the J- 
PARC E34 configuration. 

^^Here we quote the updated value from the 2010 CODATA recommended values of the 
fundamental physical constants @]. 
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component in the measurement of cja, where 21 times the data from BNL will 
be required. Equal systematic uncertainties of 70ppb are budgeted for both Ua 
and cjp, corresponding to reductions by factors of 3 and 2, respectively from what 
has been achieved. While BNLE821 improved on the CERNIII experiment in 
a revolutionary manner—primarily by the invention of direct muon injection 
into the storage ring—the ENALE989 experiment, in contrast, will introduce a 
broad suite of refinements focussed on optimizing the beam purity and rate, the 
muon storage efficiency, and modernizing the instrumentation used to measure 
both uja and ujp. 

A limiting factor at BNL was the 120 m beamline between the pion pro¬ 
duction target and the storage ring. Because the decay length of a 3.11 GeV/c 
pion is ^ 173 m, the beam injected into the storage ring contained both muons 
and a significant number of undecayed pious, the latter creating an enormous 
burst of neutrons when intercepting materials. Their subsequent capture in 
scintillator-based detectors impacted detector performance adversely. The Eer- 
milab accelerator complex will deliver pure, high-intensity muon bunches to the 
storage ring at a fill-rate frequency increase of ^ 3 compared to BNL. Proton 
batches at 8 GeV from the Booster are divided into short bunches in the Recy¬ 
cler Ring. Each bunch is then extracted and strikes a target station tuned to 
collect 3.1GeV/c 7r+ and transport them along a 270 m F0DC[3 lattice. The 
highest energy decay muons in the 7r+ ^ decay chain are captured and 

transported along the same beamline with a longitudinal polarization of ^ 97%. 
These muons (together with pious and protons) are injected into the repurposed 
p Delivery Ring (DR). There, they make several revolutions to reduce the pion 
contamination by decay and to separate protons by their velocity difference. 
A kicker in the DR extracts the pure muon bunch into a short beamline that 
terminates at the storage ring entrance. 

The centerpiece of the experiment is the 1.45 T superconducting storage 
ring ma that was re-located to Eermilab from Brookhaven in 2013. It has 
recently been reassembled in a custom temperature-controlled building having 
a firm foundation for the magnet support; both are critical for the magnetic 
field stability. Three storage ring subsystems—the superconducting inflector, 
four electric quadrupoles, and a fast kicker—determine the fraction of incoming 
muons Cstore that become stored, and their subsequent beam properties such 
as betatron oscillations. Improvements and replacements in these devices are 
aiming a factor of 2 or more increase in estore compared to BNL. 

The storage ring magnetic field will be shimmed following procedures devel¬ 
oped for E821, with small improvements owing to the need for a more highly 
uniform final field. Retooling of the pulsed NMR probes, modern 3D OPERA 
model guidance, upgrades to the in-vacuum shimming trolley, and NMR probe 
readout using waveform digitizers, represents just some of the work. The ab¬ 
solute NMR probe is the same one used in the muonium hyperfine experiment 


alternating pairs of focusing and defocusing quadrupoles magnets tuned to transport the 
pions and the forward-decay muons 
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Figure 21: Left: E821 anomalous precession data, including fit. The data is wrapped around 
every 100 fis. Right: The final, azimuthally averaged magnetic field contours; the scale is in 
ppm with respect to the average. Figures courtesy E821 collaboration. 


that established the muon-to-proton magnetic moment ratio A in Eqn.J^ and 
it will be cross calibrated with the new J-PARC experiment (see Sec. 1^ New 
absolute probes are also being developed as cross checks. 

The entire suite of detectors, electronics, calibration and data acquisition 
systems will be new and modern. These include highly segmented lead-fluoride 
Cherenkov calorimeters, where each crystal is read out by silicon photomul¬ 
tipliers uni- The signals are recorded using custom 800 MHz, 12-bit-depth 
waveform digitizers. A distributed laser-based calibration system is designed to 
maintain gain stability at the 0.04% level. Several stations of in-vacuum straw 
tracking detectors will provide beam dynamics measurements that are required 
to control various systematic uncertainties. Finally, a high-speed, GPU-based 
DAQ is prepared to process Gbytes/s of data in a deadtime free operation. 

The systematic uncertainty estimates are based on the considerable experi¬ 
ence of E821 and the targeted experimental upgrades that have been made to 
address each one of them. At the time of this review, the storage ring is being 
commissioned, with a 9-month shimming period planned for 2015. The beam- 
lines are being built, and the detector and electronics systems have undergone 
various test-beam runs at Fermilab and SLAG. Physics data taking is expected 
to begin in 2017. 

5.3.2. J-PARC muon (g-2) experiment 

While the initial 5a^ statistical precision goal of the proposed J-PARG ex¬ 
periment was also at the 100 ppb level, current estimates based on the muon 
source intensity and polarization now aim at ^ 400 ppb, similar to the BNL 
experiment |114j . The collaboration is developing a creative new method that 
will feature different systematic uncertainties compared to the storage ring ex¬ 
periment; thus, the systematic error budget is difficult to anticipate at this 
time. 

The most striking difference will be the use of a 10-times lower momentum, 
but ultra-eold^ muon beam. This choice is motivated by the desire to elim- 


50 






inate the electric field—and thus the spin-precession terms that are affected 
by it—used for vertical focusing in storage rings. To accomplish this goal, a 
muon beam must be made with a negligible transverse momentum component, 
APtIP^ ~ 10“^. This is accomplished by accelerating from rest a source of 
ultra-cold muons created by the re-ionization of muonium atoms in vacuum. 
How the overall experiment is then designed, both upstream—creating the muon 
source—and downstream—where the positrons are measured—is quite unique. 

A 3GeV proton beam strikes a graphite target, producing pions that can 
decay at rest near to the surface of the target. The 28.4MeV/c, 100% longitu¬ 
dinally polarized /i+ surface muon beam is directed to a thin target optimized 
to stop muons and form muonium atoms, M = nPe~. The target is designed 
to permit muonium to diffuse into the vacuum on the downstream surface. Ag¬ 
gressive efforts over the past few years have realized promising results in raising 
the net yield of muonium that emerges into the vacuum. A recent study at 
TRIUMF used a silica aerogel target where micro-channels in the target were 
created using a laser ablation technique HU]. This target resulted in an 8 -fold 
improvement compared to the previous measured yield using the same material. 
When the J-PARC beam rates are combined with the TRIUMF muonium yield 
measurements, an expected production of 0.2 x lOVs is found. While a factor 
of 5 lower than originally planned, it still represents a major step forward and 
a viable rate for a {g — 2) experiment. Studies will continue aimed at further 
rate improvements. 

The muonium atom ESI and its hyperfine structure measurements are de¬ 
scribed in Sec. we also direct the reader to Fig. which illustrates the 
external field dependent quantities discussed herein. Muonium formed in the 
weak magnetic field limit can be described in terms of its total angular momen¬ 
tum and associated magnetic quantum numbers (F, where the triplet and 

singlet combinations are given as: (1, l)i, (1, 0 ) 2 , (1, — 1)3 and ( 0 , 0 ) 4 , and the 
subscript i is our shorthand label for the four states. If muonium is formed 
in zero magnetic field, the relative population of the four states i = I to 4 is: 

^,0 and respectively. Allowing for a weak field, and choosing the axis of 
quantization along the incoming muon polarization direction, z, the net muon 
polarization as a function of time is given by EH] 



Here, the field strength is traditionally expressed by which is a ratio of the 
sum of the electron and muon Zeeman interactions to the muonium ground-state 
hyperfine interval, all expressed in terms of frequencies. To obtain a sample of 
at-rest, but polarized muons, the muonium atoms must be re-ionized, which 
is accomplished in the vacuum by two simultaneous laser bursts having wave¬ 
lengths Ai = 122 nm and A 2 = 355 nm. The first excites the IS to 2 P transition 
and the second ionizes the atom, leaving the free /i+ essentially at rest in the 
vacuum. Muons ionized from state I, the (1,1) triplet, retain their initial polar¬ 
ization. In contrast, those populating states 2 and 4 can make transitions at the 
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rate z /24 ~ 4.5 GHz, which is so high that the emerging muon spin will be, for all 
practical purposes, unpolarized. Consequently, the maximum net polarization 
from the muonium-at-rest source is Pz = 50%. 

The ultra-cold, “at rest,” muons will be rapidly accelerated by a linac to a 
longitudinal momentum of = 300 MeV jc. K novel feature of E34 is that the 
polarization direction can be flipped at production, prior to acceleration, by use 
of a low magnetic field. This feature might become important for anticipated 
systematic uncertainties. 

The muon beam will be directed through the top of a highly uniform 3T 
MRI-type magnet. A spiral injection path allows the muons to enter the field 
region at a steep angle, which softens such that the beam eventually orbits a 
plane perpendicular to the magnetic field, see Fig.[^ Custom fringe-fields and 
a vertical kicker are key elements in the design. Very weak magnetic focusing is 
required, but it causes a negligible perturbation to the frequency. 

Apart from the kinematic differences and uniqueness of the source, once the 
muons begin to circulate in the field, the experiment is much like the higher- 
energy storage ring designs. The muon spin processes proportionally to [g — 2) 
and the anomalous precession frequency is encoded in the modulation of rate 
vs. time of the higher-energy positron decays. 

The muon momentum and field strength values imply that the orbit radius 
is 33 cm and the cyclotron period is 7.4 ns. Decay positrons curl to the inside 
of the central orbit where vanes of silicon strip detectors are positioned. The 
radius of the reconstructed positron tracks provides the momentum (energy) 
determination with good acceptance, and the location of the detectors is tuned 
for good performance for the higher-energy range that optimizes the FOM. 
The detector system must have a stable acceptance over the measuring period 
and withstand a total initial hit rate approaching 10^ hits/s |117j . Sorting 
of hits into tracks presents a unique challenge here. However, if solvable, the 
systematics of complete track reconstruction could be lower than those inherent 
using calorimeter techniques, and in any case, they will be different. 

5.S.S. Comparison of methods 

A comparison of many of the parameters between the two new {g — 2) ex¬ 
periments is shown in Table Equation [29| can be used to evaluate the number 
of required events necessary to meet the desired statistical precision. We used 
a polarization of 50% for the J-PARC muonium beam and 97% for the decay- 
in-flight beam at Fermilab. 

5.4- Electric Dipole Moment 

A permanent electric dipole moment (EDM) of a particle or fundamental 
system violates the discrete symmetries of parity (P) and time-reversal (T). 
Because quantum field theories are CPT invariant, a T-violation observation 
leads to CP violation (CPV). Finding a new source of CPV is a major quest in 
atomic, nuclear and particle physics because of its implications in any resolution 
of the baryon-antibaryon asymmetry problem: {ub —nB^jn^^ = 6 x 10“^^ excess 
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Figure 22: The proposed setup for the J-PARC {g — 2) Experiment. Muons enter at the top 
left (green trajectory) and spiral into the highly uniform magnetic field region. Their decay 
positrons curl inward to an array of silicon tracking detectors. Figure courtesy T. Mibe. 


Table 4: Comparison of various parameters for the Fermilab and J-PARC (g — 2) Experiments 


Parameter 

Fermilab E989 

J-PARC E24 

Statistical goal 

100 ppb 

400 ppb 

Magnetic field 

1.45 T 

3.0T 

Radius 

711 cm 

33.3 cm 

Cyclotron period 

149.1 ns 

7.4 ns 

Precession frequency, uJa 

1.43 MHz 

2.96 MHz 

Lifetime, 

64.4 /is 

6.6 /is 

Typical asymmetry, A 

0.4 

0.4 

Beam polarization 

0.97 

0.50 

Events in final fit 

1.5 X 10^1 

8.1 X 10^1 
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baryons per photon in the present universe, where ub tlb = 0 is assumed at 
the Big Bang. The three Sakharov conditions required to arrive at an excess 
of baryons are: 1) at least one 5-number violating process; 2) a source of C- 
and CP-violation; and 3) interactions that occur outside of thermal equilibrium. 
The challenging task of knitting together these ingredients into a quantitatively 
complete explanation is not yet complete and the reader is referred to discus¬ 
sions of this fascinating topic; see nmns]. One important fact, is that the 
CP violation that occurs in the lone phase in the CKM mixing matrix is insuf¬ 
ficient by many orders of magnitude to explain a mechanism for baryogenesis. 
Consequently, a new source of CPV is required. Vigorous experimental efforts 
are underway at quark-flavor factories, in neutrino oscillation experiments, and 
in EDM searches using atoms, molecules and neutrons. 

We begin with a caveat. In a hadronic system, an EDM can be accom¬ 
modated within the SM owing to the CP-violating 0qcd term. The non¬ 
observation of a neutron EDM with dn < 10“^^e-cm, implies 0qcd < 10“^^. 
The SM-allowed term appears to be very finely tuned—the strong CP problem— 
if it is finite at all. The situation can be completely resolved if the physical 
axion exists; the Peccei-Quinn mechanism allows 0qcd = 0 at the price of a 
new particle [43]. An axion of the right mass is also highly motivated to fulfil 
the role of the missing dark matter particle and current searches by the ADMX 
collaboration will soon largely explore that parameter space completely. [43] 

New sources of CPV are not unexpected in popular extensions of the stan¬ 
dard model, such as supersymmetry, and two-Higgs doublet models (see Refs, 
and discussion in mS]). Current experimental limits are listed in Table [§ 
The non-observations are beginning to have severe consequences on all of the 
models. Eor example, they already severely constrain many supersymmetric 
CP-violating phases. In the simple lepton sector, the electron imbedded in a 
polar molecule is probed to very impressive limits. However, the fundamental 
source of any observed CPV in this kind of a system can have multiple inter¬ 
pretations as to its origin |12Qj . Is it singularly from the electron? The muon, 
on the other hand, is the only fundamental particle that can be directly tested 
and therefore easiest to interpret. It is also the only 2nd-generation particle 
being probed, which can have important implications in certain BSM models. 
In general, the sensitivity to new physics is expected to scale linearly owing to 
the mass term in the denominator of the dipole moment definition. Current de 
limits are then nearly 7 orders of magnitude more sensitive to new physics than 
the muon. 

However, there are BSM scenarios in which non-linear scaling occurs. Babu 
et al. mi presented an argument for scaling, a ^ 10^ fold enhance¬ 

ment, meaning a next-generation muon EDM search would be competitive to 
the electron. This analysis deserves an update based on new information from 
lepton-flavor-changing tan decay limits, the improved electron EDM, and direct 
LHC bounds, but it remains an intriguing consideration. 

More recently, Hiller et al. |122j proposed a supersymmetric model with 
CP violation from lepton flavor violation that can achieve rather large values 
for d ^—as large as 10“^^e- cm—in the extremes of the model parameter space. 
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Table 5: Selected EDM limits for the electron, Hg atom, neutron and muon. 


Type 

System 

EDM Limit (e-cm) 

Ref. 

Paramagnetic 

YbF 

de ~ 

= (- 2.4 ±5.9) X 

UMl 

Paramagnetic 

ThO 

de ~ 

= ^2-1 ±4.5) X 10-29 

[ESI 

Diamagnetic 

199jjg 

dA -- 

= (0.5 ± 1.5) X 10-29 

m\ 

Nucleon 

Neutron 

dfi - 

= (0.2 ± 1.7) X 10-29 

m\ 

Lepton 

Muon 

= 

= b0-l±0.9) X 10-1® 

m 


The practical range is constrained by the limits on the flavor-mixing decay 
BR(t ji^) < 10“^, which will be improved with Belle II running in the 
future; lower BR’s there imply smaller d^. 

The searches for SUSY at the colliders primarily focus on i?-parity con¬ 
serving models. Considering R-parity-violating (RPV) models opens up the 
parameter space considerably |123j . In such an analysis, it is observed that 
the muon is unique from the other systems being probed and limits as high 
as ^ 10cm can be imagined. While these very different examples do 
predict relatively “large” values for the relevant range sets a challenging 
experimental goal for muon enthusiasts as we discuss below. 


Experimental Considerations 

A typical atomic, molecular, or neutron EDM experiment involves a mea¬ 
surement of the difference in the spin precession frequency of a system subject 
to parallel and, alternatively, antiparallel magnetic and electric fields. In prac¬ 
tice, it is vital to work with the strongest possible electric fields, which are often 
found in the interior of atoms or molecules, where they greatly exceed labora¬ 
tory capabilities. For a relativistic muon circulating in a plane orthogonal to a 
pure dipole magnetic field, the situation is quite different. Here, the muon will 
feel a transverse induced motional electric field Em (x /3 x B. For the {g — 2) 
storage rings, where 7 = 29.3, the electric field strength is nearly 13GV/m! 

To understand the measurement concept, we rewrite Eqn. [^in the absence 
of an external electric field. While J-PARC will not use one at all, the focussing 
electric field at Fermilab has {E/c ^ B) and the 7 is selected to eliminate the 
affect on the precession of the spin. Thus, we have simply 


d^net = dJa djEDM =- 

m 




(31) 


The precession orientations for the magnetic and electric moments are or¬ 
thogonal, as illustrated in Fig. For a non-zero EDM, the precession plane 
would be tilted inward toward the center of the cyclotron orbit by the very small 
angle 

and the observed precessional frequency would be ujtot = \/^‘edm + 
key for the experimentalist is that the tilt of the plane is the signal. 


55 
















Figure 23: The net precession frequency is a vector sum of oJa caused by the anomalous 
magnet moment and ujedm caused by a possible permanent electric dipole moment. For a 
in the x direction, its spin would rota te counterclockwise in the x — y plane in the absence 
of an EDM. The ^ x B term in Eqn. |26| points along the y axis, orthogonal to uja- This tips 
the precession plane as shown. Note, \uoedm\ ^ |^a| in practice. It is exaggerated in the 
figure for clarity. 


Because the muon spin reverses every half period, the direction in which it is 
tipped also reverses. The consequence is that the observable will be an up/down 
modulation of the decay particles that is out-of-phase with the spin orientation 
by a factor of 7r/2. In practice, the most precise method of determining the tilt 
is a measurement of the average slope of the decays—upward vs downward— 
vs. time. The trajectories can be precisely determined using a set of tracking 
chambers. 

The BNL experiment found = (—0.1 ± 0.9) x 10“^^ e-cm using a lim¬ 
ited subset of the data and measured at just 1 of the 24 detector stations [7]. 
Both new {g — 2) experiments will be sensitive to an EDM at a level close to 
10“^^e-cm, a major improvement. In each case, they will rely on the up/down 
slope asymmetry using trackers. Here, the J-PARC experiment, which is an 
all-tracker detector, should have a greater overall acceptance compared to the 
FNAL experiment, which will feature at first only 3 tracker stations in the 24 
discrete detector positions. 

In the parasitic method described above, the rapid precession of the mag¬ 
netic moment reverses the upward and downward tipping of the spin owing to 
a possible EDM on every cycle. This leaves at most a very faint modulating 
signal at the frequency uJnet ~ ^a- In contrast, dedicated storage-ring measure¬ 
ments of EDMs have been proposed by Earley et al. |128j using the “frozen spin” 
technique. In this method, the uja precession is set to zero in crossed vertical 
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magnetic and radial electric fields for a proper selection of field magnitudes and 
muon momentum. If 7 ^ 7 ^ 0, the particle spin will gradually tip out of plane (up¬ 
ward or downward) processing about the radial electric field, greatly enhancing 
the signal. A compact muon EDM experiment designed on this principle has 
been suggested by Adelmann et al |129j . 

The net spin precession in Eqn. when 77 = 0, is frozen when 



(33) 


which occurs when the external radial electric field has the magnitude 



(34) 


Using the realistic set of parameters— Er = 640kV/m, = 125MeV/c, B = 

1 T—the authors of |129j predict a sensitivity of Sd^ « 7 x cm with 

1 year of running at PSI on an available muon beamline. The challenges of 
injection into this unique and compact device, {Rcyciotron = 42 cm), and the 
potential systematics that can lead to spin tipping from effects unrelated to 
an EDM are discussed in their paper, and also in Ref. |128j . While remaining 
many orders of magnitude behind the linearly-scaled de established limits, it is 
still an idea worth keeping alive. If we are to fully understand any non-zero 
EDM, we will need many probes to decouple the various interpretations from 
fundamental CPV sources. This small-scale experiment would serve well as a 
general demonstration of the storage-ring based EDM proposals that have been 
extended to focus on the deuteron and the proton, with promises in those cases 
of very impressive limits. 

6. Muonium hyperfine structure 

6.1. Experimental approaehes to muonium speetroseopy 

Muonium (/i+e“) is the electromagnetic bound state of a positive muon and 
a negative electron. It is a purely-leptonic, hydrogen-like atom that unlike either 
ordinary hydrogen or muonic hydrogen is completely free from the complications 
associated with the proton’s finite size and its electromagnetic sub-structure. 

Muons and electrons are spin-1/2 particles and consequently the muonium 
15' ground state has a hyperfine structure that comprises a spin E = 1 triplet 
state with three magnetic substates Mp = —1,0,-hi and a spin E = 0 singlet 
state with one magnetic substate Mp = 0. The interaction between the mag¬ 
netic dipole moments of the muon and the electron causes an energy splitting 
between the F = 0,1 hyperfine states of about 18 peV (4.5 GHz). As shown in 
Eig. [24] —in the presence of a static magnetic field—the Zeeman effect causes a 
further splitting of the hyperfine states, with the Breit-Rabi equation describ¬ 
ing the energy levels versus field strength (for example see Refs. |13Q[ I13l] b In 
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Figure 24: Breit-Rabi diagram of the energy levels of the muonium IS hyperfine states versus 
the field strength parameter xb = {gjfJ'B + where gj, g'^ are the electron and 

muon gyromagnetic ratios in muonium, are the electron and muon Bohr magnetons, 

Au is the hyperfine interval, and H is the field strength. The labels 1...4 denote the two 
Zeeman transitions (z /12 and 1 / 34 ) that are measured in the LAMPF and MuSEUM hyperfine 
structure experiments. 


the weak-field limit the magnetic interaction between the muon-electron mag¬ 
netic moments dominates and the aforementioned (F, Mp) are good quantum 
numbers. In the strong-field limit the magnetic interaction with the applied 
magnetic field dominates and the electron and muon spin projections (Mj, M^) 
are good quantum numbers. 

Muonium was identified in 1960 by the groups of Hughes at Columbia and 
Telegdi at Chicago The experiments detected the characteristic signal 

of muonium precession in a weak field by recording the high-energy positrons 
emitted from polarized muon stops in argon gasj^ In these circumstances the 
muonium atoms were directly formed in their IS ground state by electron cap¬ 
ture from argon atoms. Since this work, the formation of muonium has been 


such “weak field” experiments, using polarized muons and unpolarized electrons, the 
different (F, Mi?)-states are populated in proportions (1, +1) = 1/2, (1, 0) = 1/4, (1, —1) = 0, 
and (0,0) = 1/4. Such experiments thereby observe the precession of the muonium atoms 
formed in the (F, Mp) = (1, +1) hyperfine state. 
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observed in many other materials, as well as produced in vacuum through ther¬ 
mal emission from hot metal foils and fine silica powders. 

After the discovery of muonium much attention was focused on the hyperfine 
structure of its IS ground state. The first estimates of the hyperfine splitting 
were derived from measurements of the muonium polarization versus the applied 
fieldj^ Soon afterwards precision experiments involving microwave resonance 
techniques were used to induce transitions between different hyperfine states 
and thereby determine the hyperfine structure. The hyperfine transitions were 
detected through the associated muon spin-flip and the corresponding change 
in the decay positron angular distribution. 

In low-field resonance experiments, the transition frequency directly deter¬ 
mines the hyperfine interval, Au. In high-held resonance experiments, the de¬ 
termination of the hyperfine interval from the measurement of a single tran¬ 
sition frequency requires the use of the Breit-Rabi equation and knowledge of 
the muon magnetic moments. This issue motivated the development by DeVoe 
et al |133j of a so-called “double-resonance” technique involving the high-held 
measurement of two hyperfine transition frequencies using a single microwave 
cavity. The technique allows the concurrent determination of both the hyperfine 
interval and the muon magnetic moment. 

More recently, muonium hyperfine spectroscopy experiments have employed 
setups that alternate between two microwave fields to obtain the two hyperfine 
frequencies. They yield the best determinations of both the muon-to-proton 
magnetic moment ratio and the muon-to-electron mass ratio m^/mg— 

two fundamental constants of great importance to precision spectroscopy of 
muonic atoms. The ratio is also crucial to the extraction of the muon 

anomalous magnetic moment, a^, in the muon {g — 2) experiments (see Sec.j^. 

Note that the hyperfine interval is both measurable and also calculable with 
extraordinary precision (see Ref. [4] for details of theoretical calculations). With 
input of other fundamental constants—most importantly the fine structure con¬ 
stant and the Rydberg constant—the comparison between measured and cal¬ 
culated values of the hyperfine interval is considered a definitive test of QED 
theory in bound-state systems 

6.2. LAMPF hyperfine strueture experiment 

The most recent measurement of the muonium IS ground state hyperfine 
structure was conducted by Liu et al. m at LAMPF. The experiment mea¬ 
sured the frequencies of the two high-held, muon spin-flip transitions (Mj, M^) 


^®The muonium polarization is a function of the comparative sizes of the interaction energy 
of the applied field and the hyperfine splitting of the IS ground state. 

^^Although not discussed in detail here, the 1*S'-2S' interval in muonium is also measured. 
The experiment —utilizing Doppler-free, two-photon, pulse laser spectroscopy—yielded 

a 4 ppb determination of Av'is2S in good agreement with theory. Using the combination of 
the results from the hyperfine experiment and the 1S'-2*S' experiment, the authors obtained a 
verification of charge equality between muons and electrons of 1 + jq^- — ( —1.1±2.1) x 
10 - 9 . 
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Figure 25: Schematic diagram of the LAMPF HFS experiment mi indicating the muon 
beam, beam counter and profile monitor, krypton target, microwave cavity, pressure vessel 
and downstream high-energy positron detector system. Figure courtesy Liu et al. 


= (+1/2,+1/2) (+1/2,-1/2) and (-1/2,-1/2) (-1/2,+1/2), denoted 

respectively as z^i2 and z/34 in Fig. The experiment employed the double¬ 
resonance technique with a novel line-narrowing approach using a custom time- 
structured muon beam. The setup—including the muon beam, gas target, mi¬ 
crowave cavity and positron detector—is depicted in Fig. 

The experiment used a high rate, 100%-longitudinally polarized, surface 
muon beam derived from the 120 Hz repetition-rate, 650 /is pulse-period, LAMPF 
primary proton beam. An electrostatic kicker was used to produce a cycle of 
4 /is beam-on accumulation periods followed by 10 //s beam-off measurement 
periods. The beam extinction between accumulation periods was roughly 99%. 

The incident muons entered a large-bore, high-uniformity, 1.7 T magnet 
containing a copper microwave cavity filled with pure krypton gas. Muon 
stops formed polarized, ground state muonium-the (Mj, M^) = (1/2,—1/2) 
and ( — 1/2,—1/2) states—by electron capture from krypton atoms. A double¬ 
layered scintillator telescope recorded the high-energy decay positrons that were 
emitted downstream of the stopping target. A combination of plastic scintilla¬ 
tors and wire chambers were used for beam monitoring. 

The microwave cavity was designed to resonate at both the 1^12 transition 
frequency of 1897.5 MHz and the z/34 frequency of 2565.8 MHz. NMR mag- 
netometry using multiple fixed and movable probes was used to monitor the 
magnetic field. 

When precisely tuned to z/12 or z/34 the microwave field induces muon spin-flip 
transitions and thereby changes the angular distribution of decay positrons. In 
earlier measurements all positrons are detected, both those from “fast decays” 
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where the muon interacts only briefly with the microwave field and those from 
“slow decays” where the muon interacts at length with the microwave field. 
Liu et al utilized the 4 /rs beam-on, 10 //s beam-off, time structure to only 
observe the decay positrons from long-lived muonium with extended microwave 
interactions. This procedure narrowed the resonance lineshape by roughly a 
factor of three. 

The positron data was collected in super-cycles of ten beam pulses with the 
microwave cavity alternately switched between on and off and the microwave 
frequency alternately tuned to z/12 and z/34. The positron signal for deriving the 
transition frequencies was defined as 

S{u,H) = {Non/Noff-l) 

where Non {^off) represents the positron counting rates with the microwave 
field on (off) and H and u are the static field strength and the microwave 
frequency, respectively. Near the Zeeman resonances the spin-flip transitions 
caused a large increase in 

Resonance curves were collected by (i) sweeping the static magnetic field at a 
fixed microwave frequency and (ii) sweeping the microwave frequency at a fixed 
static magnetic field. The resonance curves 5'(z/, H) were then fit to lineshapes to 
extract the transition frequencies. The fitted lineshapes incorporated the muon 
stopping distribution, static field distribution, microwave power distribution, 
and positron detection efficiency. Data were collected at two gas pressures 
and extrapolated to zero pressure; the procedure accounted for a slight shift 
of the resonance frequency due to muonium-atom collisions and the resulting 
wavefunction distortions. 

The measured resonance curves determined the two Zeeman frequencies to 
precisions of 17-18 ppb (the statistical errors on the S{v^H) resonance curves 
were the dominant experimental uncertainties). Using the Breit-Rabi Eqn. and 
the values for z/12 and z/34 a hyperfine interval of 

Az/ = 4 463 302 765(53) Hz (12 ppb) 
and a magnetic moment ratio of 

= 3.183 345 13(39) (120 ppb) 

were obtained. The results were a three-fold improvement over the earlier ex¬ 
perimental work. 

As mentioned earlier, an improved value for the mass ratio m^jme can 
be obtained from the measured value for the ratio with the input of 

the muon anomalous magnetic moment and the proton-to-electron magnetic 
moment ratio jip /The precise measurement of fip/fip thus yielded a precise 
determination of the mass ratio 

mp/rrie = 206.768 277(24) (120 ppb) 

The hyperfine structure of 15' muonium thus renders the best determinations of 
both the magnetic moment ratio [ipj[ip and the lepton mass ratio mplrrie. 
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The current status of theoretical calculations of Au is given in Ref. [4]. Al¬ 
though the largest theoretical uncertainties arise from recoil correction terms, 
the overall uncertainty in the Ai/ prediction is the aforementioned experimen¬ 
tal knowledge of the mass ratio m^jme g 

The corresponding calculated and 
measured values of the hyperfine interval are in good agreement within their re¬ 
spective uncertainties of 272 Hz (61 ppb) and 53 Hz (12 ppb). This agreement 
is considered an important verification of bound-state QED calculations; the 
precision being much greater than analogous comparisons in ordinary hydrogen 
and muonic hydrogen. 

Alternatively—by equating the theoretical expression and measured value 
for the hyperfine interval Au and regarding m^/rrie as a free parameter—an 
indirect determination of the muon-to-electron mass ratio and the muon-to- 
proton magnetic moment ratio 


m^/rrie = 206.768 2843(52) [25 ppb] 

11^/lip = 3.183 345 107(84) [26 ppb] 


(35) 

(36) 


can be obtained [4]. The result for /i^//ie is important in the determination of 
the muon anomalous magnetic moment (see Sec. 5.2 and Ref. |135j b 


6.3. MuSEUM hyperfine strueture experiment 

An improved hyperfine spectroscopy experiment (MuSEUM) is under devel¬ 
opment at J-PARC |136j . The experiment will employ the same basic approach 
as the LAMPE HES experiment with the measurement of the same Zeeman 
frequencies in a high magnetic field via microwave excitation of muon spin-flip 
transitions. The experimental setup is shown in Eig. 

As mentioned the LAMPE experiment was limited by statistical uncertain¬ 
ties. The J-PARC experiment plans to increase by one order-of-magnitude the 
muon beam intensity and increase by two orders-of-magnitude the recorded de¬ 
cay positrons. To achieve these goals the experiment will combine the J-PARC 
H-line beamline designed for rates of 1 x 10^ /^^/s and a high-rate decay-positron 
detection system. 

The MuSEUM setup will incorporate a low-mass, high-rate, 2D-imaging, 
scintillating fiber hodoscope to enable pulse-by-pulse measurement of beam pro¬ 
files. The experiment will use a 1.7 T superconducting magnet with 1 ppm 
homogeneity and 50 ppb NMR magnetometry. A longer microwave cavity will 
increase the stopping fraction and reduce the gas pressure correction. A thin 
scintillator viewed by an image intensifier and CCD camera will determine the 
muon beam profile. A finely segmented, high rate, scintillator array with SiPM 
readout will provide the measurement of the downstream-going, high energy, 
decay positrons. 


^^The corrections to the hyperfine splitting from hadronic vacuum polarization and 
exchange are much smaller than the uncertainty arising from the mass ratio 171^1 me- 
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Figure 26: Schematic diagram of the J-PARC muonium experiment showing the pulsed 

muon beam, beam profile monitor, krypton gas chamber and microwave cavity, and the ab¬ 
sorber and segmented detector for high energy positron detection. Figure courtesy MuSEUM 
collaboration. 

During 2014 the MuSEUM collaboration have conducted tests of detector 
sub-systems and performed measurements of beam properties. The beginning 
of data taking is anticipated for 2015. 

7. Muonic Lamb shift 

At first blush, one might wonder why we include a discussion of the proton 
charge radius in a muon physics review. Indeed, the inclusion makes sense here 
because it is a precision muon experiment that stirred things up with a result 
now known as the “Proton Radius Puzzle” |137j . The CREMA collaboration 
at PSI reported a very precise measurement of the muonic hydrogen Lamb shift 
I/is; i.e., the 2Si/2 ^ ‘^^ 1/2 energy level difference in the /ip atom. The S energy 
levels are sensitive to the proton finite size, owing to their spherical symmetry 
and consequent overlap with the distributed charge distribution. 

The motivation for the experiment was based on the need for a more precise 
determination of the proton charge radius Vp. Along with the Rydberg constant 
Roo: h is one of the two required inputs to calculate hydrogen energy levels 
using QED, where the S-state energy is given approximately by 



(37) 


with n the usual principal quantum number. The connection to Vp can be 
obtained from the energy shift of an S-state level by 


AE = ^Tra | ^'s(O) r^. 


(38) 
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with ^s(O) the electron wavefunction at the origin. The experimental situation 
in hydrogen spectroscopy had achieved a precision great enough such that im¬ 
proved knowledge of Vp in computing the finite-size effect was limiting |137[I138] . 

In muonic-hydrogen, the Bohr radius is ^ 186 times smaller than the cor¬ 
responding one in ordinary hydrogen, which implies a greater overlap with the 
nucleus by a factor proportional to the cube of the radii, or 6.4 x 10^. Thus, 
the sensitivity to the finite-size effect is greatly enhanced. The logical sequence 
to an improved QED test is as follows: 

1. Measure the energy levels precisely in the fip system. 

2. Extract the proton charge radius Vp from the shift in S-state levels relative 
to the essentially unaffected P states. 

3. Use the improved knowledge of Vp to better compare the measured ep 
atomic energy levels to QED predictions. 

At least, that was the idea. 

However, the results of the pp measurements indicated a 4% smaller charge 
radius, with a 0.6% uncertainty, compared to what had been commonly as¬ 
sumed. The standard methods had been low-energy e — p scattering and ordi¬ 
nary hydrogen spectroscopy, where the finite size effect enters and one might 
assume QED to obtain the finite-size level shift, rather than test QED by ex¬ 
ternally knowing the finite size. The discrepancy between the muonic and 
electronic methods is significant—7 standard deviations—such that any hope 
of simply using the independent Vp extracted from muonic atoms is problem¬ 
atic. Consider the numeric values. The previous CODATA recommended 
value was Vp = 0.8775(51) fm [4], while the muon result alone gives Vp = 
0.84087(39) fm |139[ 1140] : note the 13 times smaller uncertainty in the muon 
measurement. 

Besides the statistical significance, it has been generally agreed that the 
method employed in the muonic hydrogen measurement is nearly irrefutable; 
it involves laser spectroscopy with accurately calibrated absorption lines. The 
checking and double checking of the extrapolation of Vp from AE2S-2P has re¬ 
vealed no error, or even much wiggle room of uncertainty. Perhaps the muon 
behaves differently than an electron. Perhaps it has a more complex sensitivity 
to the proton charge distribution, which would totally violate our earlier as¬ 
sertions of lepton universality. These ideas and other exotic suggestions have 
been discussed, and often dismissed or ruled out, in a vast literature that is 
summarized in Refs. |137[ 1138] . 

If there is no problem with the muon measurement, and not finding a cred¬ 
ible exotic origin for the difference, then possibly there is a problem with both 
ordinary hydrogen spectroscopy and the low-energy elastic electron-proton scat¬ 
tering measurements. These latter two methods agree with one another on r^, 
albeit with much larger error bars than the muonic Lamb shift measurement. 
Nevertheless, investigations have been raised about each method in searching 
for a resolution to the puzzle. 
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7.1. Proton radius from electron scattering and hydrogen spectroscopy 

While we refrain from a deep departure into e — p scattering formalism, the 
basic assumptions follow; more complete reviews with details, history, and un¬ 
certainty discussions exist |137[ 1138] . The relativistic electron-proton scattering 
cross section da/dQ is usually expressed in terms of a combination of the proton 
electric and magnetic form factors, Ge and Gm- These both depend separately 
on the four momentum transfer squared. One can map out the cross section 
over a series of energies and angles, and extract Ge{Q‘^) and Gm{Q‘^) using a 
so-called Rosenbluth separation This workhorse technique should serve 

well here, although one of the major findings in the JLab program has been the 
deviation of the ratio of Ge/Gm vs. at high > 3GeV^ from the Rosen¬ 
bluth method compared to a modern recoil polarization technique. Generally, 
the discrepancy is attributed to the non-inclusion of important two-photon ex¬ 
changes in the Rosenbluth method that turn out to be very important there, but 
do not impact the recoil method |142[ 1143] . We mention this only to illustrate 
that surprises can lurk in extrapolations that seem otherwise straightforward. 

The definition of the charge radius of the proton from e — p scattering is 

One must extrapolate the Ge vs. results to the unmeasurable intercept 
at zero momentum transfer. Important to the procedure is the application of 
radiative corrections, which must be under control at the sub-percent level in the 
overall error. In the most recent results from Mainz |144] . the A1 collaboration 
obtains the electron scattering (ES) proton charge radius: rp(ES) = 0.879(8) fm. 

Turning next to ordinary hydrogen spectroscopy, one recalls that only the 
S orbitals are affected by the proton finite size because they overlap with the 
charge distribution of the proton, see Eq. This means any transitions to S 
levels that are measured provide input to r^. Garlson nicely explains the problem 
of correlations in the data: the high precision value of the Rydberg constant is 
obtained from the very same atomic energy level experiments that measure the 
proton radius [138] . Percent level measurements of several transitions yield 
results that are sensitive to the finite size, but are by no means competitive to 
the precise measurements in the muon system. Gollectively the charge radius 
from hydrogen spectroscopy (HS) is: rp(HS) = 0.8758(77) fm [4]. 

7.2. Muonic Lamb shift experiment 

To study the muonic-hydrogen energy levels, one first needs to form the 
atom itself. In practice, this was done in the GREMA experiment by beginning 
with a very low-energy negative muon beam at PSI, derived from pious spiraling 
toward the center of a so-called cyclotron trap. About 30% of the pions decay 
in the 7r“ ^ /i“ + process, emitting negative muons that are further decel¬ 
erated to keV energies by multiply passing through a metalized foil at -20 kV. 
Once confined, they can leave along the axis of the trap into a toroidal mag¬ 
netic momentum filter and then onward to a target region that contains pure 
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Figure 27: The CREMA experimental layout. The low-energy muon beam enters from the 
left and is detected by the emission of electrons in a thin set of carbon foils that then excite 
scintillators. The muon passes through a E x B velocity selector to separate it from the 
electrons. It then enters an ultra-low pressure hydrogen gas target. The PMT signals indicate 
the to time and signal the laser to fire. The 5 ns long pulse is reflected in the multipass cavity 
giving an effective overlap time with the muonic atoms of about 75 ns. If the delayed 2P state 
is formed, large-area avalanche photodiodes placed near the target are positioned to record 
the characteristic 1.9 keV xrays. Figure courtesy R. Pohl. 


hydrogen gas at low pressure, see Fig. A muon entering the target will have 
passed through thin stacks of carbon foils. Low-energy electrons can be ejected, 
and subsequently detected in thin scintillators viewed by photomultiplier tubes. 
The signals so produced serve as a time marker of the arriving muon. The muon 
slowed in the process by normal dEjdx energy loss has a probability to atomi¬ 
cally capture in an excited n level—typically 14—and begin an ordinary cascade 
to the ground state. For optimized conditions of pressure, ~ 1% of the atoms 
will result in a muon cascade terminating in the 2S metastable state. This is the 
required starting point for a Lamb-shift measurement, and it was by no means 
an obvious or easy situation to prepare. For the CREMA conditions of a 1 mbar 
pressure hydrogen target, the metastable 2S state will have a lifetime of about 
1 /is before undergoing collisional de-excitation |137j . It is in that short window 
that the rest of the experiment must then work. 

Once formed, the next challenge is to induce the 2S ^ 2P level transition 
by shining a properly tuned laser on the atoms. If the laser frequency corre¬ 
sponds to AE2S-2P, the muon will transition to the appropriate 2P level, where 
it will then rapidly cascade to the IS level, emitting a characteristic 1.9 keV 
xray. These four steps are shown in the left panel of Fig. The idea of the 
experiment is to form the atom and then fire the laser 0.9 /is later, simultane¬ 
ously opening up a 75ns wide gate to observe whether a 1.9keV xray has been 
emitted. The laser frequency is tuned in discrete steps to scan the anticipated 
energy region around the expected 2S ^ 2P transition energy. The recorded 
number of xrays corresponding to the 2P ^ IS xray vs. time will have two 
features. First, a large prompt peak will be present in 99% of the cases owing to 
the atom following a normal cascade to the IS ground state. A second, delayed 
peak, at ^ 100 times lower intensity will be present if and only if the laser fre¬ 
quency is correctly tuned to the 2S ^ 2P resonance; otherwise, the second peak 
is absent. By repeating the experiment and slowly sweeping the laser frequency. 
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Figure 28: a) The four step schematic of how the experiment works. 1) The muon arrives 
in the target and captures in a high n shell. 2) Muon cascade, with ~ 1% stopping in the 
metastable 2S state. 3) A tunable laser is fired 0.9/rs after each muon arrives in the target. 4) 
The 2S —)■ 2P transition is induced if the laser frequency is on resonance, and the 2P state will 
decay rapidly to the IS ground state, emitting a xray that can be detected, serving to tag 
the process, b) Blowup of the level scheme showing the fine-structure splitting of the P levels 
and the hyperfine structure of S and P levels. The two measured and published transitions 
are indicated. 

one obtains a very precise measurement of A^ 2 S- 2 P- 

In practice, the P orbitals are split by atomic fine structure into the P 1/2 and 
P 3/2 levels, and both S and P levels are further split by the hyperfine interaction. 
The level scheme is shown on the right panel of Fig. where the Lamb shift 
is defined as shown, and the effect of the finite-size (fs) correction on the S 
shell is highlighted. The two measured transition energies, AEt (triplet) and 
AEs (singlet) can be used in linear combinations that, together with known 
QED-based corrections not dependent on r^, yield both the Lamb shift and 
the hyperfine splitting. In turn, one can not only deduce the discussed charge 
radius, but also the Zemach radius, r^, which is essentially a measure of the 
magnetic distribution inside the proton. While the extracted Vp is significantly 
more precise than other methods, rz is not. Its value of rz = 1.082(37) fm 
is compatible with other methods and its uncertainty is many times larger. 
As quoted above, we obtain here from muon spectroscopy (MS): rp{MS) = 


0.84087(39) fm 


7.3. Present and future work involving muons 

The CREMA collaboration has completed, but not yet published, additional 
Lamb-shift measurements in deuterium and helium. Three transitions in deu¬ 
terium have been accurately measured, with preliminary interpretations that 
suggest consistency with the muonic hydrogen result; however, theoretical work 
continues so it is premature to draw firm conclusions. They have also com¬ 
pleted Lamb-shift measurements on ^He and ^He systems in 2014, which are 
undergoing analysis. 
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Suppose we do not find a ready solution. Then what? A proposal by the 
MUSE collaboration |145j to help add a different data set to this discussion has 
been approved at PSI and the design of the experiment is in progress. The idea 
is to measure low-energy /i—p scattering and, parasitically, e—p scattering in a 
large-angle, open-geometry spectrometer. This novel idea would then complete 
the set of determinations of Vp from electron and muon scattering and electron 
and muon spectroscopy. If the muon is somehow different from the electron, or 
if scattering and spectroscopy methods differ, then data of this type can help 
resolve the puzzle. 

At the present time, the proposed plan aims at uncertainties roughly at the 
level of current e—p scattering measurements. Systematic uncertainties should 
be controlled by using both /i+ and beams, with their accompanying e+ 
and e~ components. Incoming tt, p and e particles in the secondary beamline 
are tagged on an event-by-event basis using time-of-ffight with respect to the 
RF structure. Because of the muon’s higher mass compared to the electron, 
radiative corrections in p—p scattering are relatively small and therefore under 
better control than for electrons, an advantage here for the muon experiment. 
Of course the difficulty of carrying out a precision scattering experiment utilizing 
a secondary decay beam introduces a complexity that will certainly challenge 
the team. 


8. Nuclear muon capture 

8.1. Basic features of muon capture 

Muon capture and beta decay are close cousins. Both processes 

1. p [Z^ A] ^ [Z — 1^ A] Up 

2. [Z,A] ^ [Z ±1,A] -h -h z/e 

involve transmutations of protons into neutrons or vice-versa through a semi- 
leptonic weak interaction with a precisely-known leptonic current. However the 
energy release in the two reactions—set by the muon mass in the muon process 
and the nuclear mass difference in the beta process—are quite different. Con¬ 
sequently, the two processes can illuminate different features of the underlying 
weak nucleonic and nuclear interactions. 

Muon capture occurs from the IS ground state of a muonic atom; such 
atoms are formed when muons are stopped in matter. In light nuclei, where 
the overlap of the muon orbital with the nuclear volume is relatively small, the 
capture rate is small compared to muon decay. In heavier nuclei, where the 
overlap is much larger, the capture rate is large compared to muon decayp^ For 
muonic hydrogen and muonic deuterium about 0.1% of muons undergo capture. 


^^The Z'^-law for muon capture ma states the capture rate is proportional to the fourth 
power of the effective charge Z of the atomic nucleus. 
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When muonic atoms are formed on non-zero spin nuclei (/ 7^ 0) the IS 
ground state is split into two distinct hyperfine states with total angular mo¬ 
menta ofF = / + l/2 and F = I — 1/2. The possibility of muon capture from 
the singlet and triplet hyperfine states in hydrogen, and doublet and quartet hy¬ 
perfine states in deuterium, is responsible for engendering capture on hydrogen 
isotopes with additional richness and additional complexity. 

The first observation of muon capture on hydrogen was reported in 1962 
by Hildebrand using a hydrogen bubble chamber. This experiment— 

together with other early experiments using bubble chambers and liquid scin¬ 
tillator detectors—were important as evidence in support of the nascent V-A 
theory of the weak interaction IMIEI]. Many other muon capture experiments 
have since been conducted. 

Muon capture on hydrogen isotopes offers a unique opportunity to determine 
elusive components of weak nuclear interactions—the induced pseudoscalar cou¬ 
pling of the proton and the two-body axial current of the deuteron. Herein we 
describe the recent progress in precision /ip and /id experiments that address 
such elementary features of weak interactions. 

8.2. Muon capture on hydrogen, jip nu 

Muon capture is generally treated as a current-current weak interaction 
where the leptonic current and the nucleonic current have the familiar par¬ 
ity violating V-A structures. The leptonic current has the simple 7^(1 —75) 
form. The nucleonic current—because of its quark constituents and their strong 
interactions—is more complicated. 

The most general form of the nucleonic V-A current is 

+9vr + 

(40) 

where 7^ are the Dirac matrices, q = Pn — Pp the momentum transfer and 
and tun are the muon and nucleon mass. The nucleonic current contains 
six “coupling constants” that are functions of the momentum transfer-squared 
. The couplings and gs are the vector, weak magnetism and induced 

scalar couplings of the hadronic vector current, V. The couplings ga^ gp, and gt 
are the axial, induced pseudoscalar and induced tensor couplings of the hadronic 
axial current, A. 

The terms involving 9m and gp are called first-class currents while the 

terms involving gs and gt are called second-class currents EHl. This distinction 
arises as first-class currents and second-class currents have opposite transfor¬ 
mation properties under G-parity—an operation that links the transmutation 
of protons into neutrons with the transmutation of neutrons into protons. Con¬ 
sequently, the second-class contributions to the leading first-class currents only 
arise through G-parity breaking effects {e.g. the u-d quark mass difference and 
the electromagnetic corrections). No experimental evidence for second-class 
currents exists (for recent discussions see Ref. |149j regarding gs and Ref. |15Qj 
regarding gt). 
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In the conserved vector current hypothesis (CVC) [21] the weak vector cur¬ 
rent and isovector electromagnetic current are the components of a conserved 
vector-isovector current. The hypothesis relates the and terms of the 
weak vector current to the charge and magnetism terms of the isovector elec¬ 
tromagnetic current. It predicts a “weak magnetism” analogous to magnetic 
effects in electromagnetic interactions and yields g^ = 1.0 and gm = 3.706 at 

— Q xhe roots of CVC were the closeness (1-2%) of the constant Gp deter¬ 
mined in muon decay and the constant Gpgv determined in beta decay. This 
observation resembles the equality of the electric charge of the electron and the 
proton. Apparently, the bare weak vector charge of the proton, like the bare 
electric charge of the proton, is protected from renormalization via emission 
and absorption of virtual pions, by a conservation law nn]. Nowadays, the 
conserved vector-isovector current is an integral part of our understanding of 
the nucleon’s quark structure. 

Concerning the axial current, the axial coupling is well-determined from 
measurements of neutron beta decay that yield ga = 1.2723 ± 0.0023 [43]. In 
contrast with g^ the value of ga is (slightly) modified by strong interactions and 
we speak of a partially conserved axial current in place of an exactly conserved 
vector current. This partially conserved axial current reflects an underlying 
approximately conserved chiral symmetry of strong interactions |152j . 

The remaining term is the induced pseudoscalar coupling gp] an interaction 
that plays a significant role in muon capture but not in beta decay. For fifty 
years the value of gp has been uncertain and the predictions for gp have been 
untested. The interest in gp stems from more than just its status as the poorly- 
known piece of the weak nucleonic current. A precise value m 

gp = 8.44 ± 0.23 (41) 

is predicted by quantum chromodynamics—a prediction that is closely con¬ 
nected to spontaneous symmetry breaking in low energy QCD and the dynam¬ 
ical origins of the hadronic masses. 

To understand the low-energy realization of chiral symmetry it is helpful to 
consider the hypothetical limit of massless u and d quarks. For rriu = = 0 

QCD possess an exact 811(2)1,xSU(2)i^ chiral symmetry, z.e., there are sepa¬ 
rate copies of isospin symmetry for the left-handed quarks and the right-handed 
quarks. This symmetry generates two conserved currents, a vector current cor¬ 
responding to the sum of left- and right-handed quark currents and a conserved 
axial current corresponding to the difference of left- and right-handed quark 
currents. At low energies this chiral symmetry is spontaneously broken through 
QCD interactions. As a result the hadrons acquire mass and the pion appears 
as the Goldstone boson of the broken symmetry. Still the underlying currents 
remain conserved currents and thereby dictate a precise relation between ga and 
9p' 

Of course, up and down quarks are not exactly massless and consequently 
chiral symmetry and axial current conservation are also not exact. This small 
explicit breaking of chiral symmetry modifies the relation between ga and gp 
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Figure 29: Schematic diagram showing the important atomic and molecular states and 

transition rates for muon stops in isotopically pure hydrogen. The fip atoms are initially 
formed in a statistical mix of triplet atoms (3/4) and singlet atoms (1/4). is the density- 

dependent triplet-to-singlet transition rate, is the density-dependent ortho-molecular 

formation rate, and Aop is the density-independent ortho-to-para transition rate. 


but still—through older current algebra techniques or newer chiral perturbation 
theory—a precise prediction for Qp results. The value of gp is thereby tied 
to our modern understanding of the strong interaction that incorporates its 
approximate chiral symmetry and partial axial current conservation as well as 
the dynamical origins of the hadronic masses (for further details see Ref. |154[ 

[US]). 

8.2.1. Muon chemistry in pure hydrogen 

Although the theoretical relation between the gp nv capture rate and the 
weak coupling constants is quite straightforward—a complication exists. The 
gp atoms that form when muons are stopped in hydrogen are small and neutral. 
Consequently, they scatter off and react with the surrounding H2 molecules, thus 
causing the F = 0,1 hyperfine populations to evolve with time. This evolution 
depends on the thermalization of the ‘hot’ gp atoms in the H2 environment 
as well as chemical reactions that form muonic molecules (the bound-states of 
a single negative muon and two hydrogen nuclei). A detailed knowledge of 
the relevant atomic and molecular processes—as shown in Fig. [29]— is therefore 
needed to extract gp from experimental data. 

The gp atom is initially formed in an excited state with a principal quantum 
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number n ^ 14. The excited atom rapidly de-excites through combinations of 
Auger emission, radiative decays and Coulomb de-excitation. On reaching the 
IS ground state the /ip-atoms have kinetic energies of typically 1 eV (for details 
see |156p and a statistical population of the hyperfine states (3/4 triplet atoms 
and 1/4 singlet atoms). 

These energetic atoms are rapidly thermalized by elastic and spin-flip col¬ 
lisions with the atomic nuclei of the surrounding H2 molecules. When their 
energies fall below the 0.18 eV /ip hyperfine splitting, the singlet-to-triplet tran¬ 
sitions are energetically forbidden and triplet-to-singlet transitions depopulate 
the higher-lying triplet state. The triplet lifetime is about 0.1 ns in liquid H2 
and about 10 ns in 10 bar H2 gas. 

At sufficient densities pp/i molecules form. Like ordinary H2 molecules, there 
exists both para /r-molecular hydrogen with a total nuclear spin / = 0 and ortho 
/i-molecular hydrogen with a total nuclear spin 1 = 1. The para molecule is the 
true ground state of the pp/i molecule. Importantly, the two molecules have 
different /ip-spin decompositions; the para-molecule being 3:1 triplet-to-singlet 
and the ortho-molecule being 1:3 triplet-to-singlet. 

The pp/i molecules are formed by collisions between /ip atoms and surround¬ 
ing H2 molecules via Auger emission /ip + iL ^ pp/i + e. The calculated rate 
of the El Auger transition forming ortho-molecules (p x 1.8 x 10^s“^ 

is much faster than the EO Auger transition forming para-molecules Xpp^^ — 
(j) X 0.75 X 10^s“^ {(j) is the H2 density normalized to the liquid H2 density 
(po = 4.25 X 10622 atoms/cm^). A recent measurement |157j of the total rate 
of molecular formation found Xp^p = p x 2.01 ±0.07 x 10^5“^ in reasonable 
agreement with theoretical predictions. 

Naively, the A/ = 0 selection rule for El transitions forbids the decay of or¬ 
tho molecules to para molecules. However, as recognized by Weinberg, through 
relativistic effects that mix ortho- and para-states the ortho pp/i molecules do 
gradually decay into para pp/i molecules. The decay rate was computed by 
Bakalov et al |158j to be Xop = 7.1 ± 1.2 x 10^s“^ (this rate is indepen¬ 
dent of density). Unfortunately, the two published measurements for Xop of 
(4.1 ± 1.4) X 10^s“^ |159j and (10.4 ± 1.4) x 10^s“^ |16Qj . are in significant 
disagreement. 

8.2.2. Experimental approaehes to (ip eapture 

The “neutron approach” to studying /ip ^ nn capture involves stopping 
muons in hydrogen and detecting the resulting 5.2 MeV capture neutrons. 
Such experiments were conducted in liquid hydrogen and gaseous hydrogen and 
achieved precisions of roughly 10% in the effective capture rate for the rele¬ 
vant F = 0,1 populations. The neutron method was, however, limited by the 
necessary determination of the neutron detection efficiencies. 

The “lifetime approach” to studying /ip ^ nn capture avoids directly detect¬ 
ing the reaction products of muon capture. Rather, it determines the capture 
rate A from the difference between the disappearance rates of the /ip atom and 
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the positive muon, z.ep^ 

where X^p = l/r^p and A^+ = l/r^+ are obtained from the measured time 
distributions of the decay electrons and positrons, respectively. The experiment 
is difficult as the decay rate is roughly 1000 times the capture rate and therefore 
the two disappearance rates are very similar—thus requiring extraordinarily 
precise /ip and /i+ lifetime measurements. The lifetime approach was pioneered 
by Bardin et al. at Saclay. 

A serious concern for both approaches was muon stops in Z > 1 surrounding 
materials. The detection of capture neutrons or decay electrons from muon 
stops in surrounding materials would alter the time distribution and distort the 
measured lifetime. Therefore the target vessel, etc., were typically constructed 
from high-Z materials so stopping muons were rapidly absorbed. Similarly, 
the small size and neutrality of /ip atoms exposes the muon to transfer to any 
Z > 1 contaminants in the H 2 gas. Again detection of capture neutrons or 
decay electrons from gas contaminants would distort the measured lifetime. 

8.2.3. MuCap experiment 

The MuCap experiment [34] was conducted at PSI. It used a custom-built, 
muon-on-demand beam to increase the sample of decay electrons while miti¬ 
gating the effects of muon pileup. It also used an active target to verify the 
stopping of muons in hydrogen and monitor the effects of gas impurities. 

The experiment was performed in H 2 gas of high chemical and isotopic purity 
at 10 bar pressure and room temperature. Under these conditions the triplet 
atoms are short-lived and the ppp molecules are rarely formed—thus preparing a 
nearly-pure sample of singlet atoms and enabling an unambiguous measurement 
of the pp singlet capture rate As. 

The experimental setup is shown in Fig. A series of incoming muon 
detectors that consisted of a plastic scintillator and a planar multiwire cham¬ 
ber together determined the arrival time and provided the pile-up protection of 
beam muons. A series of outgoing electron detectors that consisted of consecu¬ 
tive layers of cylindrical multiwire chambers and segmented plastic scintillators 
together determined the times and trajectories of the decay electron. On iden¬ 
tifying a muon the upstream electrostatic kicker was turned on and thereby the 
muon beam was turned off. 

The incoming muons were stopped in a hydrogen-filled time projection cham¬ 
ber (TPC). The TPC comprised a 5.04 liter active volume with a vertical drift 
field and a horizontal readout plane of perpendicular anode wires and cath¬ 
ode strips. The analog signals from anode wires and cathode strips were fed 
to three discriminator thresholds that triggered on: incoming muons (denoted 
EL)^ the Bragg peak of stopping muons (denoted EH)^ and the high ionization 


^"^For lip atoms the disappearance rate Xpp is the sum of the p decay rate and the p capture 
rate whereas for positive muons the disappearance rate is simply the p decay rate. 
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Figure 30: Cutaway diagram of the MuCap experiment indicating the muon counters (/iSC, 
p,PC), H 2 time projection chamber (TPC), and electron counters (ePCl, ePC2, eSC). The 
muon counters determine the muon arrival time, the electron counters determine the muon 
decay time, and the H 2 TPC validates the muon stopped H 2 gas. 
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of charged products from muon capture on gas impurities (denoted EVH). The 
discriminator hits were recorded by multihit TDCs. 

The experiment employed custom isotope separation and gas purification 
units. Isotopically pure H 2 gas was prepared from commercial, isotopically- 
pure hydrogen by repeated cycles of fractional distillation. Chemically pure 
H 2 gas was maintained by recycling the gas through a purification system that 
incorporated a cold trap and micro-porous filters. The experiment achieved a 
deuterium contamination of <10 ppb and a water contamination of about 9 ppb. 

A crystal oscillator was used for the timebase of the readout electronics. 
The collaboration was blinded to the exact frequency of the timebase during 
the data taking and the data analysis. Only after completing all the analyses 
was the frequency revealed. 

The experiment accumulated about 1.2 x 10^^ negative muon decays from 
pure H 2 gas, 0.6 x 10 ^^ positive muon decays from pure H 2 gas, as well as decay 
electrons from impurity-doped gas that permitted the investigation of muon 
transfer to gas impurities and formation rates of muonic molecules. 

The decay curves were constructed from the measured time difference {te — t^) 
between an incoming muon signal in the muon scintillator and an outgoing elec¬ 
tron signal in the electron scintillators. The TPC was used to validate that the 
muon had stopped in H 2 gas. The TPC data showed stopping muons as a trail 
of EL hits that led to several EH hits at the stop location. The algorithm 
for authenticating a stop was optimized to handle the effects of (i) the hard 
scattering of incident muons from target protons into surrounding materials, 
and (ii) the possible interference between the incoming muon ionization and the 
outgoing electron ionization in the TPC. 

The experiment also pioneered the in-situ measurement of gas impurities 
with the TPC. Muon transfer to gas impurities and subsequent capture on Z > 1 
nuclei was identified by single, delayed EVH hits at the stopping location. From 
the measured rates and time distributions of muon stops with subsequent EVH 
hits the necessary corrections due to N 2 / H 2 O impurities in the pure H 2 gas 
were then derived. 

The measured decay curves were fit to determine the muonic hydrogen life¬ 
time. In principle—due to the muon kinetics and the time evolution of the /ip 
spin-states—the theoretical time distribution is not exactly a single exponential 
decay curve. However, in practice a single exponential was a good fit to the 
time distribution and adequately determined the muon disappearance rate. 

Two corrections were necessary to extract the singlet capture rate from 
the difference — A^+ between the /i^ disappearance rates. One correction 
accounted for the small population of ppji molecules with singlet atoms. An¬ 
other correction accounted for the slight difference in the decay rate of a bound 
muon versus a free muon. After these corrections of about 18 s“^ and 12 s“^ 
respectively the final result of 

A 5 = 715.6 ± bA{stat) ± b.l{syst) 5 “^ 

was obtained Unlike many earlier experiments, the MuCap result is es¬ 

sentially free from ambiguities associated with muonic molecule formation. 
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Figure 31: The induced pseudoscalar coupling gp versus the ortho-to-para molecular tran¬ 
sition rate Xop- It shows the recent result of the MuCap experiment derived from ordinary 
muon capture in 10 bar gas and earlier results of Ref. mu for ordinary muon capture in 
liquid hydrogen and Ref. for radiative muon capture in liquid hydrogen. The data points 

indicate the experimental results iT^iTsni and the theoretical calculation of the ortho- 
to-para transition rate. The MuCap result from 10 bar H 2 gas is essentially free of ambiguities 
associated with Xop. 
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Using the latest theoretical calculations of /ip nv capture m and incor¬ 
porating radiative corrections |163j , the MuCap measurement of singlet capture 
thereby determines the coupling Pp. Using pa = 1.2701 ± 0.0025 for the axial 
coupling |164j . the MuCap result for implies a value of 


gp = 8.06 ± 0.48 ± 0.28 


for the induced pseudoscalar coupling (the uncertainties are associated with the 
MuCap measurement and the yPT calculation of the capture rate, respectively). 

The MuCap result for gp is in good agreement with the original predictions of 
current algebra and the modern predictions of chiral perturbation theory (Eqn. 
41). As shown in Fig. [sT] —the result for gp is also essentially free from the am- 
biquities associated with the limited knowledge of the ortho-to-para molecular 
transition rate that afflicted earlier experiments in liquid hydrogen. The result 
verifies our modern understanding of approximate chiral symmetry and partial 
axial current conservation in QCD. Such concepts are the foundations of our 
understanding of the origins of the hadron masses and the pion’s significance as 
the Goldstone boson of a broken symmetry. 


8.3. Muon capture on deuterium, pd nnu 

At a basic level the atomic nucleus is more than just an assembly of neutrons 
and protons. It incorporates such non-nucleonic degrees-of-freedoms as virtual 
pious and delta particles. With gy, ga, and gp all well-measured, the weak 
interaction offers a precise probe for exploring such exotic constituents of nuclear 
matter. 

Muon capture on deuterium is nature’s bridge between weak nucleonic and 
nuclear interactions. As such it parallels the role of radiative capture np dy 
on hydrogen and photo-disintegration yd ^ np of deuterium for electromagnetic 
processes. The np dy reaction provided the first unequivocal evidence for 
non-nucleonic degrees-of-freedom in electromagnetic interactions. These non- 
nucleonic effects were surprisingly large with pion currents contributing roughly 
10% of thermal neutron capture. 

The interest in exchange currents in weak interactions is more than theoret¬ 
ical. The /id ^ nnn reaction is closely related to other A = 2 weak processes 
including pp den thermonuclear fusion in stars and nd interactions in heavy- 
water neutrino detectors. It represents the only A = 2 weak interaction that is 
measurable and calculable to high precision. As such, the reaction is crucial to 
quantitatively understanding the non-nucleonic contributions to weak nuclear 
interactions and their influence on such processes as big-bang nucleosynthesis 
and stellar evolution as well as ordinary and double ;d-decay. 

The /id ^ nnn process is dominantly an allowed Gamow-Teller transition 
from the ^Si deuteron ground state to the ^Si nn continuum state. Due to the 
V-A structure of the weak interaction the deuterium capture rate from doublet 
{F = 1/2) atoms is much larger than quartet {F = 3/2) atoms. Given our 
excellent knowledge of the nucleon weak couplings and the deuteron nucleonic 
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wavefunction, the major uncertainty in jj^d capture is the poorly-known contri¬ 
bution of the two-body axial current. A precision measurement of /id capture 
can resolve this two-body current and thereby advance our theoretical under¬ 
standing of many weak nuclear processes. 

Until fairly recently the theoretical work on /id capture was based on phe¬ 
nomenological potential models of nucleon-nucleon interactions. Using phe¬ 
nomenological potentials, sophisticated calculations that incorporated detailed 
initial- and final-state nucleonic wavefunctions augmented by simplified models 
of non-nucleonic contributions were performed by Tatara et al usa, Doi et 
al. |166j and Adam et al. |167j . These calculations gave rates Kjj for doublet 
capture of typically 390-400 s~^. The calculations suggested a two-body axial 
contribution of roughly 5% that mostly originates from delta excitation via pion 
exchange between nucleons. 

Chiral effective field theory (yEFT) was first developed for systems of pions, 
then extended to systems involving a single nucleon, and eventually applied to 
few-body nuclear systems. Its development has profoundly altered the theo¬ 
retical treatment of weak interactions on few-body nuclei. yEFT established 
a rigorous, unified framework for calculations that obeys the underlying sym¬ 
metries of quantum chromodynamics while utilizing the pions and nucleons as 
low-energy degrees-of-freedom. It is based on a systematic expansion in small 
parameters—the momentum transfer, pion mass and nuclear binding energy— 
where leading-order terms are computed and higher-order terms are neglected. 
Each of the calculated terms involves a low-energy constant that must be de¬ 
termined from data. 

In yEFT, a single low energy constant determines the two-body axial current 
contribution in weak nuclear processes (this low energy constant is denoted as 
(Ir or Lia in different versions of effective theories). The /ad nnu process 
offers an unmatched opportunity for determining this constant to a precision 
comparable to the recent caculations. 

Over recent years a number of calculations of Ajj have been performed with 
increasing sophistication in the EFT framework. The two most recent calcula¬ 
tions, which consistently treat the nuclear wavefunctions and the weak opera¬ 
tors, were conducted by Marcucci et al. HIS], yielding 399T3 s ^ , and by Adam 
et al. |169j . yielding 383.8-392.4 s~^. 

8.3.1. Muon ehemistry in pure deuterium 

Just as atomic and molecular processes can complicate the interpretation 
of jip capture data, such atomic and molecular processes also complicate the 
interpretation of fid capture data. A detailed knowledge of relevant atomic and 
molecular processes—as shown in Fig. [32]— is therefore required. 

The /id atoms are formed in excited states that rapidly de-excite to the 
IS ground state by Auger emission, radiative decays and Coulomb collisions, 
thus yielding a statistical mix of “hot” doublet and quartet atoms. In many 
respects the chemical reactions of /id atoms are very similar to jap atoms. Both 
/id and /ip are tiny, neutral atoms that easily penetrate the electronic clouds 
of surrounding molecules to scatter off and react with atomic nuclei. Like /ip 
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Figure 32: Schematic diagram showing the important atomic and molecular states and 

transition rates for muon stops in isotopically pure deuterium. The fip atoms are initially 
formed in a statistical mix of quartet atoms (2/3) and doublet atoms (1/3). (l^Xqd is the 
density-dependent quartet-to-doublet transition rate and and 0A^^^ are the density- 

dependent molecular formation rates from the quartet and doublet states. Also shown is the 
muon recycling following ddfi molecule formation and /^-catalyzed fusion. 
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atoms, the jj^d atoms are thermalized by elastic and spin-flip collisions with 
surrounding nuclei. When the fid energy falls below the 0.043 eV hyperflne 
splitting, the spin-flip collisions then depopulate the quartet atoms in favor of 
doublet atoms. However, the cross sections are considerably smaller for jad d 
scattering than jap p scattering and consequently the quartet pd atoms in D 2 
are longer-lived than triplet pp atoms in H 2 (for details see |170j b 

One new feature of pd chemistry is temperature-dependent resonant forma¬ 
tion of ddp molecules—a process by which the ddp binding energy is absorbed 
by D 2 vibro-rotational modes. For example, at cryogenic temperatures, while 
ddp formation by doublet pd atoms involves a rather slow, non-resonant process, 
the ddp formation by quartet pd atoms involves a fast, resonant process. 

Another new feature of pd chemistry is muon catalyzed fusion nn]. The pos¬ 
sibility of muon catalysis of nuclear reactions was first proposed by Frank ma 
in 1947 and later considered by Gerstein, Sakharov and Zeldovich in the early 
1950s as a possible energy source. Its first observation was entirely accidental— 
Alvarev et al. [173] identifying the puzzling tracks following muon stops in bub¬ 
ble chambers as muons released following catalyzed fusion. 

This release of muons from ddp molecules is an additional dimension of pd 
chemistry. In ddp molecules the fusion reactions are 

1. dd ^ n ^He 

2. dd — y p ^H, 

the former yielding an intense source of mono-energetic neutrons in pd exper¬ 
iments. On forming a ddp molecule, the fusion reaction occurs essentially in¬ 
stantaneously. The reactions generally release muons but on occasion they will 
stick to the charged products of the fusion reaction. 

8.3.2. MuSun experiment 

The /id doublet capture in pure deuterium was previously measured using 
the lifetime technique and a liquid D 2 target yielding A^ = 470±29 s“^ |174j and 
using the neutron technique and a gaseous D 2 target yielding = 409 ±40 s“^ 
nza. These experiments were conducted more than twenty five years ago. 

The MuSun experiment [38] is using the lifetime technique to measure the pd 
doublet capture rate Ad to about 1.5% and thereby improve by roughly five-fold 
the current knowledge of two-body axial current contributions to A = 2 weak 
nuclear processes. The approach requires the preparation of a nearly-pure sam¬ 
ple of doublet atoms, a 10 ppm measurement of the pd atom lifetime, and careful 
monitoring of isotopic and chemical impurities in deuterium. The experiment 
builds on the development and the innovations in the MuCap experiment. 

As shown in Fig. the MuSun experiment is using a novel cryogenic D 2 
time projection chamber. The temperature of 34 K and pressure of 5-6 bar 
were chosen to prepare an optimal population of nearly-pure doublet atoms. At 
this temperature and pressure the quartet atoms rather quickly decay to dou¬ 
blet atoms but regeneration of quartet atoms by muon recycling following ddp 
molecule formation and dd-fusion is quite small. The time projection chamber 
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Figure 33: Cutaway diagram of the MuSun low temperature, high pressure, D 2 time projec¬ 
tion chamber. It shows the cryogenic pressure vessel, the beryllium beam entrance window 
and liquid neon cooling system, as well as the horizontal cathode plane, 6x8 segmented anode 
plane and field-shaping wires. 


also provides the stop definition for incoming muons and monitoring of /id ki¬ 
netics including both the formation of muonic molecules and the transfer to gas 
impurities. 

The MuSun setup consists of an incoming muon counter package, outgoing 
electron counter package, the cryogenic D 2 time projection chamber and a liquid 
scintillator neutron detector array. The muon and electron counter packages are 
conventional arrangements of plastic scintillators and proportional chambers 
that were originally constructed for the MuCap experiment. 

The high pressure, low temperature, D 2 TPC works as follows. Ionization is 
collected via a vertical drift field and a 6x8 segmented, horizontal anode plane, 
then readout via cryogenic pre-amplifiers and 8-bit, 25 MHz waveform digitizers. 
The TPC and associated electronics were designed for good energy resolution 
(~10 keV) to thereby enable the clean identifcation of fusion products and muon 
capture on gas contaminants. 

As in MuCap the muon and electron plastic scintillators determine the time 
interval between the incoming muon and the outgoing electron in order to con¬ 
struct the decay curve and extract the /id lifetime. A stop definition that is 
derived from the signals in the anode pads is used to validate the entries in the 
time distribution. 

Processes that involve either muon transfer to chemical impurities or molec- 
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Figure 34: Measured energy deposition of stopping negative and positive muons in the 

MuSun cryogenic TPC. The energy is defined as the summed signal from the muon stopping 
pad and its eight neighboring pads. The positive muon energy distribution shows a large 
peak that corresponds to the stopping muons and a small peak that corresponds to an in-situ 
alpha calibration source. The negative muon energy distribution additionally shows the effects 
of muon-catalyzed fusion and the additional energy deposition associated with the charged- 
particle products of the n ^He and p fusion channels. In the particular case of the p 
channel, the occurrence of one, two and three fusions following a single muon stop are clearly 
discernable. 


ular formation on isotopic impurities impose stringent limits on the possible 
chemical and isotopic contamination of the D 2 gas. Chemical impurities {e.g. 
air, water) are worrisome due to large /id^/iZ transfer rates and require concen¬ 
trations of N 2 of <1 ppb and O 2 of <3 ppb. Isotopic impurities are worrisome 
due to pdfi molecule formation and require concentrations of ordinary hydrogen 
of <10 ppm. To achieve such purities the D 2 gas was prepared in-situ by a 
custom isotope separation unit and continuously cleaned of chemical impurities 
by a custom gas recycling unit. 

The MuSun experiment is well underway at PSI. Fig. [3^ shows the measured 
energy deposition of stopping positive and negative muons in the cryogenic TPC. 
For negative muons the figure indicates the occurrence of one, two and three 
muon-catalyzed fusions are clearly detected. Production data taking was begun 
in 2013 and will continue through 2015. 


9. Summary and Outlook 

We have described many contemporary projects where muons are being used 
as probes of fundamental atomic-, nuclear- and particle-physics properties, the 
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structure of the weak interaction, and as sensitive probes of physics beyond 
the standard model. But, why muons, say, compared to electrons or tauons? 
The answer involves comparing known parameters such as mass, lifetime, and 
decay modes, as well as practical issues such as source yields, polarizability, and 
once again, lifetime. Typically, greater mass provides enhanced sensitivity to 
high-energy-scale physics. Access through radiative quantum loops generically 
scales as {mi jmheavy)‘^i where mi is the lepton mass and the subscript heavy 
might correspond to the known IT, or Z, or to an unknown new particle that 
couples to leptons. For example, given equally precise measurements of the 
lepton anomalous magnetic moments, the muon will be more sensitive than an 
electron by a factor of 43,000The tau will be 300 times better still, but it 
suffers from its fleeting lifetime of 0.29 ps, which impedes this and many other 
desired studies. 

We have seen how the relatively long muon lifetime allows for the formation 
of muonic atoms, which, with their smaller Bohr radii, can be used to sensitively 
probe nucleons and nuclei. The purely leptonic, but hydrogen-like, muonium 
atom is a probe of QED and fundamental properties, such as the muon mass 
and the magnetic moment ratio to the proton. The copious production of polar¬ 
ized muons—which are too light to decay by the strong interaction—essentially 
allows for a laboratory of weak-interaction studies. Highly polarized muons, 
and the self-analyzing nature of their decay, are essential ingredients to many 
experiments. 

Modern experiments involving muons are generally nt/i-generation efforts, 
having been designed based on earlier pioneering work. The new experiments 
excel in precision and sensitivity reach, often being complemented by equally 
important theoretical improvements. Many experiments being built now will 
have the rare characteristic of “discovery” sensitivity, and an energy-scale reach 
complementary to or beyond that of the LHC collider program. 

To summarize briefly, we first recall recent accomplishments in this field, 
and then list future directions with their planned sensitivities. 

9.1. Recent accomplishments 

• Muon lifetime: r^+ has been measured to 2.2 ps, (Ippm) by the Mu- 
Lan experiment [5] . With 2nd-order electroweak corrections now com¬ 
puted [39], the fundamental Fermi Constant is determined to be = 
1.166 378 7(6) x IQ-^ GeV“^ (0.5ppm). 

• Decay parameters: Experiments at TRIUME [56l|57| and PSI [CTI58] have 
reduced the uncertainties by up to an order of magnitude on the Michel 
parameters p, and 6, and on the transverse and longitudinal polarization 


^^The electron anomalous magnetic moment has been measured Ea 2300 times more pre¬ 
cisely compared to the muon—a spectacular achievement, which determines the fine-structure 
constant—but it does not close the gap. 
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of the electron in muon decay. Global fits [saiM] have strong implica¬ 
tions on possible deviations of the V-A structure of the weak interaction. 
Results to date confirm standard model expectations. 

• Muon anomaly: The final result of the Brookhaven {g — 2) experiment [6], 
when compared to steadily improved theoretical evaluations [m [92] of 
the muon anomaly results in a > 3cr deviation, possibly indicative of new 
physics. 

• cLFV in /i ^ ey: The MEG experiment has set the world record on a test 
of charged lepton flavor violation finding BR(/i ^ ey) < 5.7 x 10“^^ [8]. 
The limits constrain new physics involving loop processes to the 10^ TeV 
energy scale, if one assumes maximal flavor mixing. 

• Muonium HFS: The LAMPF HFS experiment measured the hyperfine 
splitting of two transitions involving four levels of the 15' muonium atom in 
the high-field limit. The results m represent an important test of bound- 
state QED and established the important ratios: la^/jXp = 3.183 345 13(39) 
(120 ppb), and rriplme = 206.768 277(24). 

• Muon capture on the proton: The MuGap experiment [34] determined 

the singlet capture rate by measuring the /i“ lifetime in 

protium gas TPG and comparing it to r^+. Including updated radiative 
corrections |163j . one obtains Qp = 8.06 ± 0.48 ± 0.28, the weak-nucleon 
pseudoscalar coupling of the proton. The result confirms a fundamental 
prediction of chiral perturbation theory and concludes a fifty year effort 
to unambiguously determine the coupling Qp. 

• Proton radius: The precise measurement |140j of two 25 2P (Lamb 
shift) transitions in muonic hydrogen determines the proton charge ra¬ 
dius, Vp = 0.84087(39) fm [140] . The results—stunningly—are 7 a smaller 
than the previous world average, which was based on e — p scattering and 
ordinary hydrogen spectroscopy. The so-called “proton radius puzzle” 
remains unsolved. 

9.2. Near-term projeets 

At the time of this review a number of approved projects are actively taking 
data or are in a construction phase. These are the ones to watch for results in 
the coming years; they include: 

• cLFV in /i ^ ey: The MEG experiment upgrade of the calorimeter, 
tracker, and other systems will improve the energy resolution and tim¬ 
ing required to achieve a sensitivity goal of 4 x 10“^^ in the next 4 years. 

• cLFV in /i ^ e conversion: GOMET at J-PARG and Mu2e at Fermilab 
will measure the coherent conversion of a muon to an electron in the field 
of a nucleus at a single event sensitivity approaching 10“^^, a bold, 4 
orders of magnitude improvement. Superconducting solenoids are used 
for particle production, transport, and final spectrometer functions. 
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• cLFV in /i ^ eee: Mu3e at PSI is in a R&D phase with an approved 
plan to reach a BR sensitivity of 10“^^. Central to their success is the 
development of ultra-thin silicon tracking detectors and, for a later phase, 
the creation of a next-generation high-intensity muon beamline. 

• Muon anomalous magnetic moment: J-PARC and Fermilab experiments [90l 
[89] are being built to reach sensitivities on Sa^ of ^ 500 ppb and 140 ppb, 
respectively. They employ radically different beam delivery and storage 
techniques, and will consequently confront different systematic errors, an 
important comparison. 

• Muon EDM: Both {g — 2) experiments can parasitically collect data sen¬ 
sitive to a muon EDM, with 1 to 2 orders of magnitude improvement 
beyond the current limit, < 10“^^e-cm. Dedicated EDM storage-ring 
plans, using a frozen-spin technique are promising, but are not yet ap¬ 
proved. 

• Muonium HES: The MuSEUM experiment will explore the hyperfine struc¬ 
ture of the muonium atom using the well-established double-resonant cav¬ 
ity technique for exciting spin-flip transitions. They aim for an order of 
magnitude increase in formed muonium atoms and two orders of magni¬ 
tude increase in recorded decays compared to the most recent LAMPE 
experiment. 

• Proton radius: CREMA will publish Lamb-shift measurements on the 
deuteron, helium-3, and helium-4 systems, important data to be compared 
to their existing hydrogen measurements. The MUSE experiment plans to 
measure low-energy g—p and e—p scattering at low Q^. These data will 
provide an important missing clue in the enduring “proton radius puzzle.” 

• Muon capture on deuterium: MuSun will complete a measurement of pd 
capture to 1.5 % precision, which will provide a clean determination of the 
low-energy constant arising in the effective-field-theory description of the 
two nucleon weak axial current. The result is relevant for fundamental as¬ 
trophysics reactions, such as pp fusion and the neutrino breakup reactions 
in the SNO experiment. 

Physicists that might identify themselves with the subfield of “Muon Physics” 
form a diverse group who happen to share a common and unique probe. Many 
of the practical issues cross traditional discipline boundaries. The word “Preci¬ 
sion” in our title evidently describes many of the results described above, but 
we naturally stretch its meaning to include ultra-high sensitivity experiments 
as well, such as those involving rare decays. We trust that the many exciting 
projects described in this review will convince the reader of the prolific record 
of this eclectic Precision Muon Physics community and of the very bright future 
that lies ahead. 
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